ASYMPTOTIC DISTRIBUTION OF QUASI-NORMAL MODES 
FOR KERR-DE SITTER BLACK HOLES 

SEMYON DYATLOV 

Abstract. We establish a Bohr-Sommerfeld type condition for quasi-normal modes of 
a slowly rotating Kerr-de Sitter black hole, providing their full asymptotic description in 
any strip of fixed width. In particular, we observe a Zeeman-like splitting of the high 
multipicity modes at a = (Schwarzschild-de Sitter), once spherical symmetry is broken. 
The numerical results presented in Appendix [B] show that the asymptotics are in fact 
accurate at very low energies and agree with the numerical results established by other 

y—{ methods in the physics literature. We also we prove that solutions of the wave equation 

can be asymptotically expanded in terms of quasi-normal modes; this confirms the validity 

p ^ of the interpretation of their real parts as frequencies of oscillations, and imaginary parts 

as decay rates of gravitational waves. 
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Quasi-normal modes (QNMs) of black holes are a topic of continued interest in theoretical 
physics: from the classical interpretation as ringdown of gravitational waves [Cj to the 
recent investigations in the context of string theory |HoHuj . The ringdownQ plays a role in 
experimental projects aimed at the detection of gravitational waves, such as LIGO |Ab) . 
O See |KoSchj for an overview of the vast physics literature on the topic and [BeCaStl IKoZhl[ 

|KoZh2t IYoUcFu] for some more recent developments. 

t— i In this paper we consider the Kerr-de Sitter model of a rotating black hole and assume 

i that the speed of rotation a is small; for a = 0, one gets the stationary Schwarzschild- 

de Sitter black hole. The de Sitter model corresponds to assuming that the cosmological 
constant A is positive, which is consistent with the current Lambda-CDM standard model 
of cosmology 

A rigorous definition of quasi-normal modes for Kerr-de Sitter black holes was given using 
the scattering resolvent in |Dyl| . In Theorem 1 below we give an asymptotic description 
of QNMs in a band of any fixed width, that is, for any bounded decay rate. The result 
confirms the heuristic analogy with the Zeeman effect: the high multiplicity modes for the 
Schwarzschild black hole split. 



1 Here is an irresistible quote of Chandrasekhar ". . . we may expect that any intial perturbation will, 
during its last stages, decay in a manner characteristic of the black hole itself and independent of the 
cause. In other words, we may expect that during these last stages, the black hole emits gravitational 
waves with frequencies and rates of damping that are characteristic of the black hole itself, in the manner 
of a bell sounding its last dying notes." 
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Theorem 2 confirms the standard interpretation of QNMs as complex frequencies of 
exponentially decaying gravitational waves; namely, we show that the solutions of the 
scalar linear wave equation in the Kerr-de Sitter background can be expanded in terms of 
QNMs. 

In the mathematics literature quasi-normal modes of black holes were studied by Bache- 
lot, Motet-Bachelot, and Pravica |Bal[ \Ba2\ IBaMoBal iPr] using the methods of scattering 
theory. QNMs of Schwarzschild-de Sitter metric were then investigated by Sa Barreto- 



Zworski |SaBaZwj . resulting in the lattice of pseudopoles given by (0.3) below. For this 
case, Bony-Hafner |BoHa] established polynomial cutoff resolvent estimates and a res- 
onance expansion, Melrose-Sa Barreto-Vasy [MeSaBaVa] obtained exponential decay for 
solutions to the wave equation up to the event horizons, and Dafermos-Rodnianski |DaRolj 
used physical space methods to obtain decay of linear waves better than any power of t. 

Quasi-normal modes for Kerr-de Sitter were rigorously defined in |Dyl| and exponential 
decay beyond event horizons was proved in Dy2| . Vasy |Va] has recently obtained a 



microlocal description of the scattering resolvent and in particular recovered the results 
of |Dyl[ |Dy2| on meromorphy of the resolvent and exponential decay; see |Val Appendix] 
for how his work relates to |Dyl| . The crucial component for obtaining exponential decay 
was the work of Wunsch-Zworski |WuZwj on resolvent estimates for normally hyperbolic 
trapping. 

We add that there have been many papers on decay of linear waves for Schwarzschild 
and Kerr black hole s — se e [XnBll IBIStl rDaRo2l IDaRo3l IDoSchSoll IDoSchSo2l IFiKaSm Yal 
IFiKaSmYaErr} \Ta\ ITaTol IToj and references given there. In that case the cosmological 
constant is (unlike in the de Sitter case, where it is positive), and the methods of scattering 
theory are harder to apply because of an asymptotically Euclidean infinity. 

Theorem 1. Fix the mass M of the black hole and the cosmological constant A. (See 



Section 1.1 for details.) Then there exists a constant a > such that for \a\ < a® and each 
vq, there exist constants C w , such that the set of quasi-normal modes oj satisfying 

Reu > C w , Imw > — u (0.1) 

coincides modulo 0(|w| _oc ) with the set of pseudopoles 

u = F(m, l,k), m,l,keZ, < m < C m , \k\ < I. (0.2) 

(Since the set of QNMs is symmetric with respect to the imaginary axis, one also gets an 
asymptotic description for Re a; negative. Also, by |Dyl[ Theorem 4], all QNMs lie in the 
lower half-plane.) Here T is a complex valued classical symbo^ of order 1 in the (I, k) 



2 As in |Dyl| , the indices w,m, ... next to constants, symbols, operators, and functions do not imply 
differentiation. 

3 Here 'symbol' means a microlocal symbol as in for example |Tay[ Section 8.1]. For the proofs, however, 



we will mostly use semiclassical symbols, as defined in Section 2.1 



ASYMPTOTIC DISTRIBUTION OF QNMS FOR KERR-DE SITTER 



3 



variables, defined and smooth in the cone {m G [0, C m ], \k\ < 1} C M. 3 . The principal 
symbol J-q of T is real-valued and independent of m; moreover, 



7 = V \Jm^ + " ^ + V2)] + 0{l ~ X) ^ a = °' ( °' 3) 



(%F )(m, ±k, k) = (2 + 2 ^° 2)Q + 0(a 2 ). (0.4) 



The pseudopoles (0.2) can be computed numerically; we have implemented this com- 
putation in a special case I — \k\ — 0(1) and compared the pseudopoles with the QNMs 
computed by the authors of [B eCaStj . The results are described in Appendix [Bj One 
should note that the quantization condition of [Sa BaZw] was stated up to 0(/ _1 ) error, 
while Theorem [l] has error 0(l~°°); we demonstrate numerically that increasing the order 
of the quantization condition leads to a substantially better approximation. 



Another difference between (0.2) and the quantization condition of [SaBaZwJ is the extra 



parameter k, resulting from the lack of spherical symmetry of the problem. In fact, for 



a = each pole in (0.3) has multiplicity 21 + 1; for a ^ this pole splits into 21 + 1 distinct 
QNMs, each corresponding to its own value of k, the angular momentum with respect to the 
axis of rotation. (The resulting QNMs do not coincide for small values of a, as illustrated 



by (0.4)). In the physics literature this is considered an analogue of the Zeeman effect. 



Since the proof of Theorem [T] only uses micro local analysis away from the event horizons, 



it implies estimates on the cutoff resolvent polynomial in cu (Proposition 1.3). Combining 
these with the detailed analysis away from the trapped set (and in particular near the 
event horizons) by Vasy [Vaj . we obtain estimates on the resolvent on the whole space 



(Proposition 1.2). These in turn allow a contour deformation argument leading to an 
expansion of waves in terms of quasinormal modes. Such expansions have a long tradition 
in scattering theory going back to Lax-Phillips and Vainberg — see |TaZw] for the strongly 
trapping case and for references. 

For Schwarzschild-de Sitter black holes a full expansion involving infinite sums over 
quasinormal modes was obtained in [B oHa] (see also |ChrZw] for simpler expansions in- 
volving infinite sums over resonances). The next theorem presents an expansion of waves 
for Kerr-de Sitter black holes in the same style as the Bony-Hafner expansion: 

Theorem 2. Under the assumptions of Theorem^ take uq > such that for some e > 0, 
every QNM u has \ Imu + z/ | > e. (Such u exists and can be chosen arbitrarily large, as 



the imaginary parts of QNMs lie within 0(\a\ + I ) of those in (0.3).,) Then for s large 



enough depending on u , there exists a constant C such that every solution u to the Cauchy 
problem on the Kerr-de Sitter space 

U g u = 0, u\ t , =Q = / 6 H S (X S ), <9 t *u| t * =0 = fx e H S -\X S ), (0.5) 
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where Xg = (r_ — 5, r + + 5) x § 2 the space slice, t* is the time variable, and 5 > is a 
small constant (see Section 1.1 for details), satisfies for t* > 0, 

H s-i). (0.6) 
(0.7) 



\u(t*)-n vo (foJx)(t*)\\Hi(x s ) < Ce-"* (WUhs + \\h\\ 



Here 



E e ~ lt * Q E (ty^Afoji)-, 

the outer sum is over QNMs Q, Jq is the algebraic multiplicity ofu as a pole of the scattering 
resolvent, and Hqj are finite rank operators mapping H S (X$) © H s ~ l (X$) — > C°°(Xs). 
Moreover, for \u)\ large enough (that is, for all but a finite number of QNMs in the considered 
strip), Jq = 1, Uq q has rank one, and 



in. 



< C\uj 



N-s 



Here N is a constant depending on u , but not on s; therefore, the series (0.7) converges 
in H 1 for s > N + 2. 



The proofs start with the Teukolsky separation of variables used in Dyl , which reduces 



our problem to obtaining quantization conditions and resolvent estimates for certain radial 
and angular operators (Propositions |1.5 | a nd 1.6). These conditions are stated and used 
to obtain Theorems [T] and [2] in Section [T| Also, at the end of Section L2 we present the 
separation argument in the simpler special case a = 0, for convenience of the reader. 

In the spherically symmetric case a = 0, the angular problem is the eigenvalue problem 
for the Laplace-Beltrami operator on the round sphere. For a ^ 0, the angular operator 
Pg is not selfadjoint; however, in the semiclassical scaling it is an operator of real principal 
type with completely integrable Hamiltonian flow. We can then use some of the methods 
of |HiSj| to obtain a microlocal normal form for h 2 Pg\ since our perturbation is O(h), we 
are able to avoid using analyticity of the coefficients of Pg. The quantization condition we 
get is global, similarly to |VuNg| . The proof is contained in Section [3] it uses various tools 
from semiclassical analysis described in Section [2] 

To complete the proof of the angular quantization condition, we need to extract informa- 
tion about the joint spectrum of h 2 Pg and hD^ from the microlocal normal form; for that, 
we formulate a Grushin problem for several commuting operators. The problem that needs 
to be overcome here is that existence of joint spectrum is only guaranteed by exact commu- 
tation of the original operators, while semiclassical methods always give 0(h°°) errors. This 
complication does not appear in |HiSj[ |HiSjVuNg] as they study the spectrum of a single 



operator, nor in earlier works |Chl VuNg on joint spectrum of differential operators, as 



they use spectral theory of selfadjoint operators. Since this part of the construction can be 
formulated independently of the rest, we describe Grushin problems for several operators 
in an abstract setting in Appendix [Aj 

The radial problem is equivalent to one-dimensional semiclassical potential scattering. 
The principal part of the potential is real-valued and has a unique quadratic maximum; 
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the proof of the quantization condition follows the methods developed in jCdVPat IRal SJT 



In |CdVPaj . the microlocal behavior of the principal symbol near a hyperbolic critical point 
is studied in detail; however, only self-adjoint operators are considered and the phenomenon 
that gives rise to resonances in our case does not appear. The latter phenomenon is studied 
in |Ra] and |Sjl ; our radial quantization condition, proved in Section |4j can be viewed as a 
consequence of |Ral Theorems 2 and 4]. However, we do not compute the scattering matrix, 
which simplifies the calculations; we also avoid using analyticity of the potential near its 
maximum and formulate the quantization condition by means of real microlocal analysis 
instead of the action integral in the complex plane. As in |Raj . we use analyticity of the 
potential near infinity and the exact WKB method to relate the microlocal approximate 
solutions to the outgoing condition at infinity; however, the construction is somewhat 
simplified compared to |Ra[ Sections 2 and 3] using the special form of the potential. 

It would be interesting to see whether our statements still hold if one perturbs the metric, 
or if one drops the assumption of smallness of a. Near the event horizons, we rely on |Dyl, 
Section 6], which uses a perturbation argument (thus smallness of a) and analyticity of 
the metric near the event horizons. Same applies to Section 4.2| of the present paper; the 
exact WKB construction there requires analyticity and Proposition 4J2 uses that the values 
u± defined in (4.9) are nonzero, which might not be true for large a. However, it is very 



possible that the construction of the scattering resolvent of [VaJ can be used instead. The 
methods of |Vaj are stable under rather general perturbations, see |Val Section 2.7], and 
apply in particular to Kerr-de Sitter black holes with a satisfying |Va| (6.12)]. 

A more serious problem is the fact that Theorem [T] is a quantization condition, and thus 
is expected to hold only when the geodesic flow is completely integrable, at least on the 



trapped set. For large a, the separation of variables of Section |1.2| is still valid, and it is 
conceivable that the global structure of the angular integrable system in Section 3J2 and 
of the radial barrier-top Schrodinger operator in Section 4.1| would be preserved, yielding 
Theorem [T] in this case. Even then, the proof of Theorem 2] no longer applies as it relies on 
having gaps between the imaginary parts of resonances, which might disappear for large a. 

However, a generic smooth perturbation of the metric supported near the trapped set 
will destroy complete integrability and thus any hope of obtaining Theorem [TJ One way of 
dealing with this is to impose the condition that the geodesic flow is completely integrable 
on the trapped set. In principle, the global analysis of |VuNg| together with the methods for 
handling 0(h) nonself adjoint perturbations developed in Section [3] and Appendix |A| should 
provide the quantization condition in the direction of the trapped set, while the barrier- 
top resonance analysis of Section 4^3 should handle the transversal directions. However, 
without separation of variables one might need to merge these methods and construct a 
normal form at the trapped set which is not presented here. 

Another possibility is to try to establish Theorem [2] without a quantization condition, 
perhaps under the (stable under perturbations) assumption that the trapped set is normally 
hyperbolic as in jWuZwj . However, this will require to rethink the contour deformation 
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argument, as it is not clear which contour to deform to when there is no stratification of 
resonances by depth, corresponding to the parameter m in Theorem [TJ 

1. Proofs of Theorems [T] and [2] 

1.1. Kerr— de Sitter metric. First of all, we define Kerr-de Sitter metric and briefly 
review how solutions of the wave equation are related to the scattering resolvent; see 
also |Dyl[ Section 1] and |Vat Section 6]. The metric is given by 



dr 2 d9 2 \ A e sin 2 # 
A7 + A^V " (l + a)V' 



+ — — ■ — r (dt — a sin 2 9 dip)' 



A r 

;i +^)v 

Here 9 G [0, it] and <p G R/27rZ are the spherical coordinates on § 2 and r,t take values in 
M; M Q is the mass of the black hole, A is the cosmological constant, and a is the angular 
momentum; 

Ar 2 ~ 

y 

Aa 2 



A r = (r 2 + a 2 ) (l - - 2M r, A e = 1 + a cos 2 9, 

2 2,2 2 a 

p = r + a cos 9, a 



3 

The metric in the (t, r, 9, ip) coordinates is defined for A, r > 0; we assume that this happens 
on an open interval r G (r_ , r+) , where r± are two of the roots of the fourth order polynomial 
equation A r (r) = 0. The metric becomes singular at r = r±; however, this apparent 
singularity goes away if we consider the following version of the Kerr-star coordinates 
(see |DaRo2t Section 5.1] and |TaToj ) : 

t* = t-F t (r), <p* = <p-F ip (r), (1.1) 

with the functions F t , F v blowing up like c± log |r — r±\ as r approaches r±. One can choose 
F t , F^, so that the metric continues smoothly across the surfaces {r = r±}, called event 
horizons, to 

M 5 = R t x X s , X 5 = (r_ - 5, r + + 5) x § 2 , 

with 5 > is a small constant. Moreover, the surfaces {t* = const} are spacelike, while 
the surfaces {r = const} are timelike for r G (r_,r + ), spacelike for r G" [r_,r + ], and null 
for r G {r_,r + }. See |Dyl[ Section 1], |Dy2[ Section 1.1], or |Va[ Section 6.4] for more 



information on how to construct F t , F^ with these properties. 

Let O g be the d'Alembert-Beltrami operator of the Kerr-de Sitter metric. Take / G 
H s ~ 1 (Ms) for some s > 1, and furthermore assume that / is supported in {0 < t* < 1}. 
Then, since the boundary of Ms is spacelike and every positive time oriented vector at 
dMs points outside of Ms, by the theory of hyperbolic equations (see for example |DaRo2[ 
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Proposition 3.1.1] or |Tay[ Sections 2.8 and 7.7]) there exists unique solution u G Hf oc (Ms) 
to the problem 

U g u = f, suppu C {t* > 0}. (1.2) 

We will henceforth consider the problem (1.2); the Cauchy problem ( 0.5[ ) can be reduced 
to (1.2) as follows. Assume that u solves (0.5) with some f G H s (Xs), f\ G H s ~ 1 (Xs). 
Take a function x £ C°°(R) such that suppx C {t* > 0} and supp(l - x) C {t* < 1}; then 
x(t*)u solves ( 1.2[ ) with / = yjw supported in {0 < t* < 1} and the fP -1 norm of / is 
controlled by |/o||h* + ||/i||ff— 

Since the metric is stationary, there exists a constant C e such that every solution u 
to (1.2) grows slower than e^ 8-1 ^* ; see |Dyl[ Proposition 1.1]. Therefore, the Fourier- 
Laplace transform 



ulu) 



u(t*)df 



is well-defined and holomorphic in {Imcu > C e }. Here both u(t*) and u(u) are functions 
on Xg. Moreover, if f(oo) is the Fourier-Laplace transform of /, then 



P g (u)u(u) 



p 2 f{u), Imw > C e 



1.3) 



where P g (oo) is the stationary d'Alembert-Beltrami operator, obtained by replacing D t * 
with —u in p 2 O g . (The p 2 factor will prove useful in the next subsection.) Finally, since / 
is supported in {0 < t* < 1}, the function f(u>) is holomorphic in the entire C, and 



u} s /MHjjj-i^ (x,) <C\\f\\ H s- HMs) 



for Imcu bounded by a fixed constant. Here Hf , /i > 0, is the semiclassical Sobolev 



space, consisting of the same functions as iP -1 , but with norm || (hD) s ~ l f\\tf instead of 

iiprviu- 

If -P g (w) was, say, an elliptic operator, then the equation (1.3) would have many solutions; 
however, because of the degeneracies occuring at the event horizons, the requirement that 
u G H s acts as a boundary condition. This situation was examined in detail in |Vaj ; the 
following proposition follows from |Val Theorem 1.2 and Lemma 3.1] (see |Dyl[ Proposi- 
tion 1.2] for the cutoff version): 

Proposition 1.1. Fix uq > 0. Then for s large enough depending on uq, there exists a 
family of operators (called the scattering resolvent) 

R(u) : H'-^Xs) ->• H S (X S ), Imw > -u , 



meromorphic with poles of finite rank and such that for u solving (1.2), we have 

R(u)f(u), Imoo > C e . 



U{tU) 



(1.5) 
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Note that even though we originally defined the left-hand side of (1.5) for Imu > C e , 
the right-hand side of this equation makes sense in a wider region Imu > — u , and in fact 
in the entire complex plane if / is smooth. The idea now is to use Fourier inversion formula 

u(t*) = ^- I e- iwt *R{u)f{u)duj (1.6) 

and deform the contour of integration to {Imw = — z/ } to get exponential decay via the 
e -iu>t f ac t or _ We pick up residues from the poles of R(u) when deforming the contour; 
therefore, one defines quasi-normal modes as the poles of R(uj). 

Our ability to deform the contour and estimate the resulting integral depends on having 
polynomial resolvent estimates. To formulate these, let us give the technical 

Definition 1.1. Let h > be a parameter and 7Z(u; h) : Hi — » H.2, u £ U(h) C C, be a 
meromorphic family of operators, with Hj Hilbert spaces. Let also Q(h) C U(K) be open 
and 2(h) C C be a finite subset; we allow elements of Z(h) to have multiplicities. We 
say that the poles of 7Z in Q(h) are simple with a polynomial resolvent estimate and given 
modulo 0(h°°) by Z{h), if for h small enough, there exist maps Q and U from Z(h) to C 
and the algebra of bounded operators "Hi — > H2, respectively, such that: 

• for each Q' e Z(h), uj = Q(cu') is a pole of TZ, \u — Q'\ = 0(h°°), and II(u;') is a 
rank one operator; 

• there exists a constant N such that ||n(cD')||^ 1 _ s .^ 2 = 0(h~ N ) for each Q' e Z(h) 
and, moreover, 



^ uj - Q(u') 

Q'ez(h) v ' 



In particular, every pole of TZ in Q(h) lies in the image of Q. 

The quantization condition and resolvent estimate that we need to prove Theorems [T] 
and |2] are contained in 

Proposition 1.2. Fix u > and let h > be a parameter. Then for a small enough 
(independently ofu ), the poles of R(u) in the region 

I Ima;| < u , h' 1 < I Reu\ < 2h~ l , (1.7) 

are simple with a polynomial resolvent estimate and given modulo 0(h°°) by 

uj = h^ 1 J rw (m, hL hk\ h), m.l.k G Z, 

V ; (1.8) 

< m < C m , Cf 1 <hl< C u \k\ <l. 

Here C m and G\ are some constants and T u {m, I, k; h) is a classical symbol: 

F*(m, I, ~k;h)~Yl hjjP j( m > I ~ k )- 



ASYMPTOTIC DISTRIBUTION OF QNMS FOR KERR-DE SITTER 

The principal symbol J-~q is real-valued and independent of m; moreover, 

6\/61V1q 

(2 + 9AM 2 )a 



1 _ OA TWf 2 ~ 

" ±(l + h/2-ih{m + 1/2)) + 0{h 2 ) for a 



3y/3M 
(d~ k J^)(m,±~k,~k) 



0. 



27M 2 



0(a 2 



Finally, if we consider R(u) as a family of operato rs be tween the semiclassical Sobolev 
spaces H^ 1 — > H^, then the constant N in Definition 1.1 is independent of s. 



Theorem [T] follows from the here almost immediately. Indeed, since R(u) is independent 
of h, each is homogeneous in (/, k) variables of degree 1 — j; we can then extend this 
function homogeneously to the cone |A;| < I and define the (non-semiclassical) symbol 

F(m,l,k) ~ ^2j^(m,l,k). 

Note that F{m,l,k) = h- 1 J zu (m,hl,hk;h) + 0(h°°) whenever Cf 1 < hi < Q. We can 
then cover the region (0.1 ) for large with the regions (1.7) for a sequence of small values 
of h to see that QNMs in (0.1) are given by (0.2) modulo 0(|w| _0 °). 

Now, we prove Theorem [2] Let u be a solution to (1.2), with / e H s ~ l and s large 
enough. We claim that one can deform the contour in (1.6) to get 

1 



u(t* 



i Res a 

Im£>-i'o 



—iut* 



R(u)f(u)] 



2tt 



-iuit 



R(u)f(u) doj. 



1.9) 



Imt 



-I'D 



The series in (1.9) is over QNMs cD; all but a finite number of them in the region {Imw > 
— vq\ are equal to Q{u') for some Q' given by (0.2 ) and the residue in this case is e~ %ut ' Ii{u')f{Q) 
Here Q and II are taken from Definition 1.1. Now, by (1.4), we have 

||n(a/)/(w)||^ < cr<o?>^— 



\\R(u)f(u)\\ m <C(u) 



N-i 



H 3 - 1 1 

HB-i, Imw = — Vq, 



for some constant N independent of s; therefore, for s large enough, the series in (1.9) con- 
verges in H 1 and the H 1 norm of the integral in (1.9) can be estimated by Ce~ uot ||/ \h s - 1 -, 
thus proving Theorem [2] 

To prove (1.9), take small h > 0. There are 0(h~ 2 ) QNMs in the region (1.7); therefore, 
by pigeonhole principle we can find uo(h) G [/i _1 ,2/i _1 ] such that there are no QNMs 
h 2 -c\ose to the segments 

l±{h) = {Rea; = ±uo(h), —u < Imu < C e }. 

Then \\R(u)f(u)\\ H i = 0(h s ~ N ~ 4 ) on 7±(/i); we can now apply the residue theorem to the 
rectangle formed from 7±(/i) and segments of the lines {Imw = C e },{Imw = —^0}, and 
then let h —> 0. 
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1.2. Sep aration of variables. First of all, using |Vat (A. 2), (A. 3)], we reduce Proposi- 
tion |l.2| to the following^ 



Proposition 1.3. Take 5 > and put 

K s = (r_ + S,r + - 5) x § 2 , R g (u) = l Kt R{u)l Ks : L 2 {K S ) H 2 {K & ). 

Then for a small enough^ and fixed Vq, the poles of R g (u) in the region (1.7) are simple 
with a polynomial resolvent estimate L 2 — > L 2 and given modulo 0(h°°) by (1.8). 



Furthermore, by |Val Proposition A.l] the family of operators R g (u) 



coincides with the 



one constructed in |Dyl[ Theorem 2], if the functions Ft,F v in (1.1) are chosen so that 
(t, if) = (t*, </?*) in M t x Kg- We now review how the construction of R g (u) in |Dyl| works 



and reduce Proposition [L3] to two separate spectral problems in the radial and the angular 
variables. For the convenience of reader, we include the simpler separation of variables 
procedure for the case a = at the end of this section. 

First of all, the operator P g {ui) is invariant under the rotation ip i— )■ ip + s; therefore, the 
spaces T>' k = Ker(D (/ , — k) of functions of angular momentum k G Z are invariant under 
both Pgiuj) and R g {oj). In |Dyl| , we construct R g {oj) by piecing together the restrictions 
R g (u, k) = R g (uj)\x>i for all k. Then, Proposition 1.3| follows from 



Proposition 1.4. Under the assumptions of Proposition 1.3 , there exists a constant 
such that for each k G Z, 



norm 



(1) if h\k\ > Ck, then R g (u,k) has no poles in the region ( 1.7[ ) and its L 2 — > L 2 
is 0(\k\~ 2 ); (This is a reformulation of |Dyl[ Proposition 3.3] J 

(2) if h\k\ < Ck, then the poles of R g {uj,k) in the region (1.7) are simple with a poly- 
nomial resolvent estimate L 2 — > L 2 and given modulo 0(h°°) by (1.8), with this 
particular value of k. 



Now, we recall from |Dyl[ Section 1] that the restriction of P 9 (uj) to V k has the for™^' 

P r (u, k) + Pq(cj)\d> , where 



P r (u, k) = D r (A r D r ) - ^ + . a ^\ (r 2 + a 2 )u - ak) 2 , 



A., 



Po(") = —^D e {A e smeD 6 
sm(J 



t . h (ausm 2 9-D^) 2 
Aq sm 9 



1.10) 



are differential operators in r and (9, <p), respectively. Then R g (u, k) is constructed in |Dyl[ 
Proof of Theorem 1] using a certain contour integral |Dyl[ (2.1)] and the radial and angular 



4 One could also try to apply the results of |DaVa] here, but we use the slightly simpler construction 
of |Val Appendix], exploiting the fact that we have information on the exact cutoff resolvent. 
5 The smallness of a is implied in all following statements. 

6 The operator P 9 {uj) of |Dyl| differs from our operator by the conjugation done in |Va| Appendix]; 
however, the two coincide in Kg. 
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resolvents 

R r (uj,X,k) : L 2 omp (r_,r+) -> i?£ c (r_,r+), 

J2fl(w, A) : L 2 (§ 2 ) -> if 2 (S 2 ), A G C; 

i? r is a certain right inverse to P r (u,k) + A, while Rg is the inverse to P$(uS) — A; we 
write Rg(uj,X,k) = Rg(cj, . Recall that both R r and Rg are meromorphic families of 
operators, as defined in |Dyl[ Definition 2.1]; in particular, for a fixed value of uj, these 
families are meromorphic in A with poles of finite rank. By definition of R g (uj, k), a number 
u 6 C is a pole of this operator if and only if there exists A G C such that (uj, A, k) is a 
pole of both R r and Rq. 

Now, for small h > we put 

uj = hReuj, v = Imu, X = h 2 KeX, fi = h\mX, k = hk; (1-H) 



the assumptions of Proposition 1.4 2) imply that 1 < uj < 2, |P| < z/o, and < C^. 
Moreover, |Dyl[ Proposition 3.4] suggests that under these assumptions, all values of A for 
which (lu, A, k) is a pole of both R r and Rg have to satisfy |A|, \£l\ < C\, for some constant 
C x . 

We are now ready to state the quantization conditions and resolvent estimates for R r 
and Rg; the former is proved in Section [4] and the latter, in Section [3] 

Proposition 1.5 (Radial lemma). Let C\ be a fixed constant and put K r = (r_ +5, r + — 5). 

Then the poles of l^ r .R r (co>, A, k)l Kr as a function of X, in the region 

1 < u < 2, \v\ < i/ , < C fc , |A|, |/2| < Ca, (1-12) 

are simple with polynomial resolvent estimate L 2 — > L 2 (in the sense of Definition 1.1) and 
given modulo 0(h°°) by 

X + ihjl — J^im, 6j, v 1 k; h), m G Z, < m < C m , (1-13) 

for some constant C m . The principal part J-^ of the classical symbol J 71 " is real-valued, 
independent of m and v, and 

3\/3Mo 1 2 



ih(m + 1/2) + " Y "^^° -Ju + ihu) 
y/l - 9AM 2 

~. 27M 2 ^ 2 6a~ku 2 



+ 0(h 2 ) for a = 0, 



In particular, for co,k satisfying (1.12), every pole X satisfies X > e for some constant 
e > 0. 

Proposition 1.6 (Angular lemma). Let Cg be a fixed constant. Then the poles of Rg(u, A, k) 
as a function of X in the region 

1 < uj < 2, \v>\< u , \k\ < C fe , Cg 1 <X<Cg, < Cg, (1-14) 
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are simple with polynomial resolvent estimate L 2 — > L? and given modulo 0{h°°) by 

~X + ihft = J*(hl,Cb,i>,k\h), leZ, max(|jfe|, Cf 1 ) < hi < C h (1.15) 

for some constant C\. The principal part T§ of the classical symbol T e is real-valued, 
independent of v , and 

J d = l(J+ h) + 0(h°°) for a = 0. 

Moreover, J 7 Q(±k,Cj,k) = (1 + a) 2 (k — au) 2 , dfj^^k, u, k) = ±2k + 0(a 2 ), and conse- 
quently, dfrJ-'oi^kjUi, k) = —2au + 0(a 2 ). 



Combining Propositions 1.5 and 1.6 with the results of |Dyl| , we get 



Proof of Proposition 1.4 We let J-^m, /, k; h) be the solution u + ihv to the equation 

J-" r (m, u), z>, k; h) = T 6 (I , ui , v, k; h). (1-16) 

We can see that this equation has unique solution by writing T r — T e = T' + ihT" and 
examining the principal parts of the real-valued symbols J 7 ', J 7 " for a = 0. 

The idea now is to construct an admissible contour in the sense of |Dyl[ Definition 2.3]; 
e.g. a contour that separates the sets of poles (in the variable A) of R r and Rg from each 
other; then [Dyl], (2.1)] provides a formula for R g (u, k), which can be used to get a resolvent 
estimate. We will use the method of proof of |Dyl[ Proposition 3.4]. Take the contour 7 
introduced there, for l 2 = Cxh^ 1 , li — L — C\h~ 2 , and C\ some large constant. Then we 
know that all angular poles are to the right of 7 (in Tz)- Moreover, the only radial poles 
to the right of 7 lie in the domain {| Im A| < I2, |A| < L} and they are contained in the 
set {Ag, . . . , Xc m } f° r some constant C m , where A^(co>, k) is the radial pole corresponding 
to h- 2 T r (m, u, v, k; h). In particular, those radial poles are contained in 

U x = {C, 1 <\<C e , < C e }, 

for some constant Cg. 

Assume that u is not a pole of R g (u,k); then we can consider the admissible contour 
composed of 7 and the circles j m , < m < C m , enclosing AJ n (w, k), but none of the 
other poles of R r or Rg. Using the meromorphic decomposition of R r at \ r m and letting its 
principal part be W m /(\ — A^), we get 

R g (uj,k) = y^nr(u),k)®R4(u),\ r m (u,k),k) + — I Rr(u,\,k) ® Rg(u),\,k)d\. (1.17) 
^— ' Ztti J 

Here we only include the poles \ T m lying to the right of 7; one might need to change l\ in 
the definition of 7 a little bit in case some \ r m comes close to 7. The integral in (1.17) is 



holomorphic and bounded polynomially in h, by the bounds for R r given by Proposition L5 
together with the estimates in the proof of |Dyl Proposition 3.4]. 

Now, the poles of Rg in U\ are given by (1.15); let Af (w, k) be the pole corresponding to 
h^ 2 J rd (hl, u, u, k; h) and U® /(A— Af ) be the principal part of the corresponding meromorphic 
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Figure 1 . The contour 7 and interaction between radial poles (denoted by 
dots) and angular poles (denoted by asterisks). 



decomposition. Then the resolvent estimates on Rq given by Proposition 1.6 together 



with (1.17) imply 



Hg(U, k) - Hol(u) + > — — — — — — . 



1.1* 



Here Hoi (a;) is a family of operators holomorphic in u and bounded polynomially in h. 
Moreover, 

\ r m (u, k) - Af (w, fc) = h^iJ^im, uj, v, k; h) - J* (hi, Co, v, k; h)) + 0{h°°); (1.19) 



therefore, the equation AJ ra (w, fe) — Af(w, fc) = is an 0(h°°) perturbation of (1.16) and it 
has a unique solution u m> i(k), which is 0(h°°) close to Yr x T^{yn, hi, k; h). Finally, in the 



region (1.7) we can write by (1.18) 



ml 



with Hol(a;) as above and Ii m ^{k) being the product of a coefficient polynomially bounded 
in h with (W m ® Hf)(uj mt i(k), k); this finishes the proof. □ 
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Finally, let us present the simplified separation of variables for the case a = 0, namely 
the Schwarzschild-de Sitter metric: 

Ar, 2 r 2 



g = ^Ldt 2 - —dr 2 - r 2 (d6 2 + sin 2 9 dtp 2 ), 
Ar 2 



A r = r 2 (l-^-J -2M Q r. 

Note that d6 2 + sin 2 6 dip 2 is just the round metric on the unit sphere. The metric decouples 
without the need to take Fourier series in ip; the stationary d'Alembert-Beltrami operator 
has the form P g {uS) = P r (u>) + Pe, where 



i 



P r (u) = D r (A r D r ) - —u 



2 



A 



r 



1 D 2 
Pe = —MsiTieDe) + : 



sin# uy u ' sin 2 9 
Here Pq is the Laplace-Beltrami operator on the round sphere; it is self-adjoint (thus 
no need for the contour integral construction of |Dyl[ Section 2]) and is known to have 
eigenvalues 1(1 + 1), where I > 0. Each such eigenvalue has multiplicity 2Z + 1, corresponding 
to the values —I, . . . , I of the (^-angular momentum k. The angular Lemma 1.6 follows 
immediately (We nevertheless give a more microlocal explanation in this case at the end 



of Section 3.1 ) One can now decompose L 2 into an orthogonal sum of the eigenspaces of 
Pe; on the space V\ corresponding to the eigenvalue A = 1(1 + 1), we have 

P g (u)\ Vx = P r (co) + A. 



Therefore, the only problem is to show the radial Lemma L5 in this case, which is in fact 
no simpler than the general case. (Note that we take a different path here than [SaBaZwJ 
and |BoHaj . using only real microlocal analysis near the trapped set, which immediately 
gives polynomial resolvent bounds.) 



2. Preliminaries 



2.1. Pseudodifferential operators and microlocalization. First of all, we review the 
classes of semiclassical pseudodifferential operators on manifolds and introduce notation 
used for these classes; see |EvZwj Sections 9.3 and 14.2] or |DiSj | for more information. 

For k G K, we consider the symbol class S' fc (lR n ) consisting of functions a(x, £; h) smooth 
in (x, £) e M. 2n and satisfying the following growth conditions: for each compact set K C W 1 
and each pair of multiindices a, (3, there exists a constant C a $K such that 

\d^a(x,C; h)\ < C a p K {£) k -W, x G K, £ e R n , h > 0. 

If we treat M 2n as the cotangent bundle to MP, then the class S h is invariant under changes of 
variables; this makes it possible, given a manifold M, to define the class S k (M) of symbols 
depending on (x,£) G T*M. 
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If k G R and ctj(x, £) G S k ■'(M), j = 0, 1, . . . , is a sequence of symbols, then there exists 
the asymptotic sum 

i>o 

i.e., a symbol a(x, f ; h) G S fc (M) such that for every J = 0, 1, ... , 



a(x, 



£;/0- Va 3 (x,0eh J S k - J . 

0<j<J 



The asymptotic sum a is unique modulo the class h°°S 00 of symbols all of whose derivatives 
decay faster than h N (£) for each iV on any compact set in x. If a is given by an asymptotic 



sum of the form (2.1), then we call it a classical symbol and write a G S^(M). We call 
ao(x, £) the principal part of the symbol a(x,£; /i). 

Let \l/ fc (M) be the algebra of (properly supported) semiclassical pseudodifferential op- 
erators on M with symbols in S k (M). If H™ loc (M), m G K, consists of functions lo- 
cally lying in the semiclassical Sobolev space, then every element of \I/ (M) is continuous 
H™ loc (M) -)■ H™^ k (M) with every operator seminorm being 0(1) as /t — >• 0. Let *^(M) 
be the algebra of operators with symbols in S^(M) and $ c i(M) be the union of for all 
k. Next, let the operator class /i°°\l/ _00 (M) correspond to the symbol class h co S~ 00 (M)] it 
can be characterized as follows: A G h OD ty~°°(M) if and only if for each N, A is continuous 
H^ loc (M) — > Hfr loc (M), with every operator seminorm being 0(h N ). The full symbol of 
an element of ^/ k (M) cannot be recovered as a function on T*M; however, if A G \1/^(M), 
then the principal symbol of A is an invariantly defined function on the cotangent bundle. 
If M is an open subset of M. n , then we can define the full symbol of a pseudodifferential 
operator modulo h^S' 00 ; we will always use Weyl quantization. 

We now introduce microlocalization; see also |EvZwl Section 8.4] and |SjZw[ Section 3]. 



Define U C T*M to be conic at infinity, if there exists a conic set V C T*M such that 
the symmetric difference of U and V is bounded when restricted to every compact subset 
of M. For a G S k (M) and U C T*M open and conic at infinity, we say that a is rapidly 
decaying on U, if for every V C U closed in T*M, conic at infinity, and with compact 
projection onto M, every derivative of a decays on V faster than h (0~ N f° r every N. 
We say that A G ^ k (M) vanishes microlocally on U if its full symbol (in any coordinate 
system) is rapidly decaying on U. If A, B G \I/ fc (M), then we say that A = B microlocally 
on U, if A — B vanishes microlocally on U. 

For A G \l/ fc (M), we define the semiclassical wavefront set WF/,(A) C T*M as follows: 
^ WFh(A) if and only if A vanishes microlocally on some neighborhood of (x, £). 
The set WF h (A) is closed; however, it need not be conic at infinity. Next, we say that A 
is compactly microlocalized, if there exists a compact set K C T*M such that A vanishes 
microlocally on T*M \ K. We denote by \E' comp (M) the set of compactly microlocalized 
operators. Here are some properties of microlocalization: 
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If A vanishes microlocally on U\ and U2, then it vanishes microlocally on U\ U 1]%. 
The set of pseudodifferential operators vanishing microlocally on some open and 
conic at infinity U C T*M is a two-sided ideal; so is the set of operators with 
wavefront set contained in some closed V C T*M. In particular, WF ^(AB) C 
WF h (A)nWF h (B). 

A vanishes microlocally on the whole T*M if and only if it lies in ft 00 ^ -00 . 
If A vanishes microlocally on U, then WF h(A) D U — 0; the converse is true if U is 
bounded. However Q A does not necessarily vanish microlocally on the complement 
of WFh(A); for example, the operator A = e~ x l h lies in \1/ and has an empty 
wavefront set, yet it does not lie in ft°°\I/ _00 . 

The set v|/ comp forms a two-sided ideal and it lies in ty~ N for every N. 
Each A G \|> comp vanishes microlocally on the complement of WF^(A). 



• Let A G #ci( M ) and let its symbol in some coordinate system have the form (|2.1l); 
introduce 

V = supp Oj. 

j>o 

Then A vanishes microlocally on some open set U if and only if U (1 V = 0; A is 
compactly microlocalized if is bounded, and WF h(A) is the closure of V. 

We now consider microlocally defined operators. Let U C T*M be open. A local pseudo- 
differential operator A on U is, by definition, a map 

B^[A-B]e ^ comp (M)//i 00 ^- 00 (M), 5 G ^ comp (M), WF h (B) C C/, 

such that: 

• WF h ([A.S]) C WF h (S). 

• If B 1} B 2 G ^ com P(M) and WF^Bj) C f/, then [A ■ (B 1 + S 2 )] = [A ■ B x ] + [A ■ B 2 }. 

• If C G ^ fc (M), then [A • 5]C = [A • (SC)]. 

We denote by \& loc (?7) the set of all local operators on U. Note that a local operator is 
only defined modulo an /i°°v|/ -00 remainder. If A G \& fe (M), then the corresponding local 
operator A is given by [A ■ B] = AB mod h 00 ^^ 00 ; we say that A represents A. For 
M = W 1 and U C T*M, there is a one-to-one correspondence between local operators 
and their full symbols modulo h°°; the symbols of local operators are functions a(x,£;h) 
smooth in G U all of whose derivatives are uniformly bounded in h on compact 

subsets of U. In fact, for a symbol a(x, £; ft,), the corresponding local operator is defined 
by [A ■ B] = (a#b) w (x,hD x ), where b(x,£;h) is the full symbol of B G \I/ com P; since 6 is 
compactly supported inside U modulo 0(h°°) and a is defined on Z7, we can define the 
symbol product a#6 uniquely modulo 0(h°°). In particular, a classical local operator 



A G \I/|.° C (M) is uniquely determined by the terms of the decomposition (2.1) of its full 



7 This issue can be avoided if we consider WFh(A) as a subset of the fiber compactified cotangent bundle 



T*M, as in |Val Section 2.1]. Then an operator is compactly microlocalized if and only if its wavefront set 
does not intersect the fiber infinity. 
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symbol. Note that U is not required to be conic at infinity, and we do not impose any 
conditions on the growth of a as £ — > oo. 

Local operators form a sheaf of algebras; that is, one can multiply local operators defined 
on the same set, restrict a local operator to a smaller set, and reconstruct a local operator 
from its restrictions to members of some finite open covering of U. This makes it possible 
to describe any local operator A G ty loc (U) on a manifold using its full symbols in various 
coordinate charts. For A G \I/ loc ([/), one can define its wavefront set WF h{A) as follows: 
(^0,60) ^ WF/ l (A) if and only if the full symbol of A is 0(h°°) in some neighborhood of 
(xq,£q). If A G fy k represents A G \I/ loc , then WF h (A) = WF h (A); in general, wavefront 
sets of local operators obey WF h {A + B) C WF ft (A)UWF ft (B) and WF h (AB) C WF h (A)n 
WF h (B). 

Finally, we study microlocalization of arbitrary operators. Let Mi and M 2 be two mani- 
folds. An /i-dependent family of (properly supported) operators A(h) : C°°(Mi) — > V'(M 2 ) 
is called tempered, or polynomially bounded, if for every compact K\ C Mi, there exist 
iV and C such that for any u G C™{K\), \\A{h)u\\ H -N < Ch~ N \\u\\ H N. Note that the 

composition of a tempered operator with an element of is still tempered. We can also 
treat distributions on M2 as operators from a singleton to M2. 

For a tempered family we define its wavefront set WFft(A) C T*(M\ x M2) as 

follows: (x, £; 77) ^ WFh(A), if and only if there exist neighborhoods Ui(x, £) and U2(y, rf) 
such that for every Bj G \[' comp (M ? ) with WF ft ( J B i ) C t/^, we have B 2 A(h)B x G /i 00 ^" 00 . We 
say that -Ai = A 2 microlocally in some open and bounded U C T*M, if WF/ l (A 1 — A 2 )nf/ = 
0. Also, is said to be compactly microlocalized, if there exist Cj G \l/ comp (Mj) such 

that A(h) — C 2 A(h)Ci G /i 00 ^!/ -00 . In this case, all operator norms ||A(/i)|| „n 1 n 2 are 

h ^ h 

equivalent modulo 0(h°°); if any of these norms is 0(h r ) for some constant r, we write 
||A(/i)|| = 0(h r ). Here are some properties: 

• If A is compactly microlocalized, then WFh(A) is compact. The converse, however, 
need not be true. 

• If A G \l/ fc (M), then the two definitions of compact microlocalization of A (via its 
symbol and as a tempered family of operators) agree; the wavefront set of A as a 
tempered family of operators is just {(x, | (x, £) G WF/ l (A)}. 

• If v4i, v4 2 are two tempered operators and at least one of them is either compactly 
microlocalized or pseudodifferential, then the product A 2 Ai is a tempered operator, 
and 

WF/ l (v4 2 v4i) c WF A (A0 o WF fc (A 2 ) 
= {(z,£;*,0 I 3(2/,??) : (*,£;</, 77) G WF,(A), (y^;^C) e WF ft (A 2 )}. 
Moreover, if both Ai,A 2 are compactly microlocalized, so is A 2 A\. 



Let us quote the following microlocalization fact for oscillatory integrals, which is the start- 
ing point for the construction of semiclassical Fourier integral operators used in Section 2.3 
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Proposition 2.1. Assume that M is a manifold, U C Mxl m is open, ip(x,9) is a smooth 
real-valued function on U , with x G M and 9 G W 71 , and a(x, 6) G Cq°(U). Then the family 
of distributions 

u(x) = [ e Mx ' e)/h a(x,9)d9 



is compactly microlocalized and 

WF^(u) C {(x, d x ip(x, 9)) | (x, 9) G suppa, dg(p(x,9) = 0}. 

2.2. Ellipticity and formal functional calculus. Assume that U C T*M is open and 
A G \l/|.° c ([/). We say that A is (semiclassically) elliptic on U if its principal symbol does 
not vanish on U. Under this condition, there exists unique operator A~ l G \fr^p(U) such 
that A~ l A, AA^ 1 = I as local operators. The next proposition provides the form of the 
symbol of A -1 ; it is based on the standard parametrix construction: 

Proposition 2.2. Fix a coordinate system on M. Assume that A G V^(U) is elliptic 
and has the full symbol a ~ a + ha\ + h 2 a2 + . . . . Then A~ x has the full symbol b ~ 
bo + hbi + h 2 b 2 + • • • , where each bj is a linear combination with constant coefficients of the 
terms of the form 

M 
m=l 

Here the (multi) indices a m ,(3 m ,l m satisfy the condition 

M M M 

^ |^m| ^ \fim\ j ^ lm- 
m=l m=l m=l 

Furthermore, we can assume that \a m \ + \/3 m \ + l m > for all m. 



Proof. We call (2.2) an expression of type (M a , Mp, L), where M a} Mp,L are the sums of 
\oc m \i \Pm\, an d lm, respectively. If / is an expression of type (M a ,Mp,L), then we can 
prove by induction that d^d^f is an expression of type (M a + \at\, Mg + L). Now, we 
write the equation aj^b = 1; the principal term gives bo = a^ 1 , and the next terms give that 
each bj is the sum of expressions of type (j — L,j — L, L), by induction and the formula for 
the symbol product aj^-b. □ 

Let A G yP^; we say that it is elliptic on an open conic at infinity U C T*M in the class 
Yjyfc ^ or m icrolocally elliptic), if its principal symbol ao satisfies \ao(x, £)| > (£) k /C(K) for 
(x, £) in any given closed conic at infinity K C U with compact projection onto M, and 
some constant C(K) depending on K. In this case, the full symbol of A~ l satisfies the 
decay conditions of the class ^>~ k in U. In particular, if A G ^ is elliptic in the class ^ h 
everywhere, then we can define A -1 G ^"j^ for h small enough. 

We now construct functional calculus of local real principal pseudodifferential operators. 
For this, we use holomorphic functional calculus [DuSchl Section 7.3]; another approach 
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would be via almost analytic continuation DiSj , Chapter 8]. First, assume that A G 



$™ mp (M) has compactly supported Schwartz kernel. In particular, the principal symbol 
a of A is compactly supported; let K C C be the image of a . Let f(z) be holomorphic in 
a neighborhood Q of K, and let 7 C Q be a contour such that K lies inside of 7. For each 
h, the operator A is bounded I? — >■ L 2 ; for ft, small enough, its spectrum lies inside of 7. 
Then we can define the operator f(A) by the formula 

f(A) = 7^— / f(z)(z - A)' 1 dz. 
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For z G 7, the operator 2 — A is elliptic in the class \I/ ; therefore, (z — A)" 1 G ^(M). It 



follows that /(A) G \1/°(M). By Proposition 2.2, the full symbol of f(A) (in any coordinate 
system) is the asymptotic sum 

00 2j 

E^E^^^^w^o. (2-3) 

j=0 M=0 

Here a ~ a + /ia 1 + . . . is the full symbol of A; bjM are the functions resulting from applying 
certain nonlinear differential operators to do, cl\, ■ ■ ■ ■ 

Now, assume that U C T*M is open and A G \I/|,° C (?7) has real-valued principal symbol a . 



Then the formula (2.3 ) can be used to define an operator f[A] G ^° C (C/) for any / G C°°(R). 
Note that the principal symbol of (z — A)^ 1 is (z — ao) _1 ; therefore, the principal symbol 
of /[A] is / o ao- The constructed operation posesses the following properties of functional 
calculus: 

Proposition 2.3. Assume that U C T*M is open, A G ^(t/), and f>9 e C°°(K)- T/ien: 
i. WF/i(/[i]) C aQ 1 (supp/) ; where ao is i/ie principal symbol of A. 
2- Iff(t) = YjuJ)!' ^ a polynomial, then f[A] = f(A), where f(A) = Ej/i^'- 

3. (f + g)[A) = f[A)+g[A) and (fg)[A] = f[A]g[A]. 

4. If Be ^ C {U) and [A, B] = 0, then [f[A},B] = 0. 

The identities in parts 2 — 4 are equalities of local operators; in particular, they include 
the ft 00 !]/ -00 error. In fact, the operator f[A] is only defined uniquely modulo h 00 ^^ 00 . 



Proof. 1. Follows immediately from (2.3) 



2. Take an open set V compactly contained in U; then there exists A G \I/™ mp (M) such 
that A = A microlocally on V. Since / is entire and A is compactly microlocalized, we 
can define f(A) by means of holomorphic functional calculus; it is be a pseudodifferen- 
tial operator representing f[A]. Now, f(A) = f(A) microlocally on V by properties of 



multiplication of pseudodifferential operators and f[A] = f[A] microlocally on V by (2.3); 
therefore, f(A) = f[A] microlocally on V. Since V was arbitrary, we have f(A) = f[A] 
microlocally on the whole U. 
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3. We only prove the second statement. It suffices to show that for every coordinate 
system on M, the full symbols of and /[A] (7 [A] are equal. However, the terms in 
the asymptotic decomposition of the full symbol of (fg) [A] — f[A]g[A] at (x,£) only depend 
on the derivatives of the full symbol of A at (x, £) and the derivatives of / and g at a (x, £). 
Therefore, it suffices to consider the case when / and g are polynomials. In this case, we 
can use the previous part of the proposition and the fact that f(A)g(A) = (fg)(A). 

4. This is proven similarly to the previous part, using the fact that [A, B] = yields 
[f(A),B] = for every polynomial /. □ 

Finally, under certain conditions on the growth of / and the symbol of A at infinity, f[A] 
is a globally defined operator: 

Proposition 2.4. Assume that A e ^(M), with k > 0, and that A is elliptic in the class 
outside of a compact subset ofT*M. Also, assume that f G C°°(R) is a symbol of order 
s, in the sense that for each I, there exists a constant C\ such that 

|/ (0 (f)|<C,<t)-', teR. 

Then f[A] is represented by an operator in ty s c k (M). 



Proof. We use (2.3); by Proposition 2.2, the symbol bjM lies in S kM K Since / is a symbol 



of order s, /' m ms a symbol of order s — M. Then, since a E S k is elliptic outside of a 



compact set and k > 0, we have f( M > o a® £ S k( - S M \ It follows that each term in (2.3) lies 



in S sk J ; therefore, this asymptotic sum gives an element of \E^j f. □ 

2.3. Quantizing canonical transformations. Assume that Mi and M 2 are two man- 
ifolds of the same dimension. Recall that the symplectic form Uj on T*Mj is given by 
Uj = do-j, where cr^ = t;dx is the canonical 1-form. We let Kj C T*Mj be compact and 
assume that $ : T*M\ — > T*Mi is a symplectomorphism defined in a neighborhood of K\ 
and such that = K^. Then the form erf — $*<rf is closed; we say that $ is an exact 

symplectomorphism if this form is exact. Define the classical action over a closed curve in 
T*Mj as the integral of o^j over this curve; then $ is exact if and only if for each closed 
curve 7 in the domain of $, the classical action over 7 is equal to the classical action over 
$07. We can quantize exact symplectomorphisms as follows: 

Proposition 2.5. Assume that $ is an exact symplectomorphism. Then there exist h- 
dependent families of operators 

B 1 : V\M X ) -> C °°(M 2 ), B 2 : V'{M 2 ) CS°(Mi) 

such that: 

1. Each Bj is compactly microlocalized and has operator norm 0(1); moreover, WF/ l (5 1 ) 
is contained in the graph 0/$ and WFf,(5 2 ) is contained in the graph o/$ _1 . 

2. The operators B\B 2 and B 2 B\ are equal to the identity microlocally near K 2 x K 2 and 
K\ x K\, respectively. 
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3. For each P G \& c i(Mi), there exists Q G \l/ c i(M 2 ) that is intertwined with P via B\ 
and B 2 : 

B X P = QB U PB 2 = B 2 Q 

microlocally near K\ x K 2 and K 2 x K\, respectively. Similarly, for each Q G ty c \(M 2 ) there 
exists P G ty c \(Mi) intertwined with it. Finally, if P and Q are intertwined via B\ and B 2 
and p and q are their principal symbols, then p = q o $ near K\ . 

If the properties 1-3 hold, we say that the pair (Bi, B 2 ) quantizes the canonical transfor- 
mation <£> near K\ x K 2 . 

Proof. We take B\,B 2 to be semiclassical Fourier integral operators associated with $ and 
respectively; their symbols are taken compactly supported and elliptic in a neigh- 
borhood of K\ x K 2 . The existence of globally defined elliptic symbols follows from the 
exactness of $; the rest follows from calculus of Fourier integral operators. See |GuiSt[ 



Chapter 8] or VuNg Chapter 2] for more details. □ 



Note that the operators B\ and B 2 quantizing a given canonical transformation are not 
unique. In fact, if Xj G ^° mp (Mj) are elliptic near Kj and Yj G ^/ c c ° mp (Mj) are their 
inverses near Kj, then (X 2 BiXi,YiB 2 Y 2 ) also quantizes $; moreover, P is intertwined 
with Q via the new pair of operators if and only if X\PY\ is intertwined with Y 2 QX 2 via 
(B l7 B 2 ). 

We now study micro local properties of Schrodinger propagators. Take A G \l/™ mp (M) 
with compactly supported Schwartz kernel and let ao be its principal symbol; we assume 
that a is real- valued. In this case the Hamiltonian flow exp(tH ao ), t G M, is a family 
of symplectomorphisms defined on the whole T*M; it is the identity outside of suppa . 
Moreover, exp(tH ao ) is exact; indeed, if V — H ao , then by Cartan's formula 

Cycr s = d{iy(y s ) + iyd(a s ) = d{i v a s — a ) 

is exact. Therefore, 

d t exp(tV)*o- s = exp(tV)*C v a s 

is exact and exp(tV)*a s — a s is exact for all t. 

For each t, define the operator exp(itA/h) as the solution to the Schrodinger equation 

hD t exp(itA/h) = Aexp(itA/h) = exp(itA/h)A 

in the algebra of bounded operators on L 2 (M), with the initial condition exp(iOA/h) = 1. 
Such a family exists since A is a bounded operator on L 2 (M) for all h. Here are some of 
its properties (see also |EvZwl Chapter 10]): 

Proposition 2.6. 1. The operator exp(itA/h) — I is compactly microlocalized and has 
operator norm 0(1). 

2. If A,B G \l/™ mp have real-valued principal symbols and [A, B] = 0(h°°), then 
[exp(itA/h),B] = 0(h°°), exp(it(A + B)/h) = exp(itA/h) exp(itB/h) + 0(h°°). 
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(We do not specify the functional spaces as the estimated families of operators are compactly 
microlocalized, so all Sobolev norms are equivalent.) In particular, if B = 0{h°°), then the 
propagators of A and A + B are the same modulo 0(h°°). 

3. Let P G * c i and take 

P t = exp(itA/h)P exp(-itA/h). 

Then P t is pseudodifferential and its full symbol depends smoothly on t. The principal 
symbol of P t is p o exp(tH ao ), where po is the principal symbol of P; moreover, 

WF h (P t ) = exp(-tH ao )(WF h (P)). 

4- Let K C T*M be a compact set invariant under the Hamiltonian flow of a . If X G 
^comp ^ e g Ua if i(j, en tity microlocally near K , then the pair (X exp(—itA/h), X exp(itA/h)) 
quantizes the canonical transformation exp(tH ao ) near K x K. Moreover, if P,Q e ty c ° mp 
are intertwined via these two operators, then Q = exp(—itA/h)Pexp(itA/h) microlocally 
near K . 

5. Assume that V is a compactly supported vector field on M , and let exp(tV) : M — >■ M 
be the corresponding flow, defined for all t ; denote by exp(tV)* the pull-back operator, acting 
on functions on M. Let K C T*M be compact and invariant under the flow of V , and 
X G \I/™ mp have real-valued principal symbol and be equal to the identity microlocally near 
K; consider {hV/i)X G ^° mp . Then for each t, 

exp(it(hV/i)X/h) = exp(tV)* 

microlocally near K x K . 

The statements above are true locally uniformly in t. 

Proof. 1. First, take u G L 2 (M); then, since the principal symbol of A is real- valued, we 
have || A — A*\\ L 2^. L 2 = 0(h) and thus 

D t \\ exp {itA/h)u\\ 2 L2 = h^dA - A*) exp(itA/h)u, exp(itA/h)u) L 2 = 0(\\ exp(itA/h)u\\ 2 L2 ); 
therefore, exp(itA/h) is tempered: 

||exp(^//i)|| L2 ^ L2 = 0(e c; l*l). 
The rest follows from the identity 

it i f f i 
exp(itA/h) = I + —A + -A / (t - s) exp(isA/h) ds ■ -A. 

h h Jo h 

2. We have 

D t (exp(itA/h)B exp(-itAfh)) = h^ 1 exp(itA/h)[A, B] exp(-itAfh) = 0(h°°); 
this proves the first identity. The second one is proved in a similar fashion: 
D t (exp(-it(A + B)/h) exp(itA/h) exp(itB/h)) = 0(h°°). 
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3. We construct a family P t of classical pseudodifferential operators, each equal to P 
micro locally outside of a compact set, solving the initial- value problem 

D t P t = h- 1 [A,P t \ + 0(h 00 ), P = P + O(h°°). 

For that, we can write a countable system of equations on the components of the full symbol 
of P t . In particular, if p(t) is the principal symbol of P t , we get 

dtp(t) = {a ,p(t)} = H ao p(t); 

it follows that p(t) = po o exp(tH ao ). Similarly we can recover the wavefront set of P t from 
that of P. Now, 

dt(exp(-itA/h)P t exp(itA/h)) = 0(h°°); 

therefore, P t = P t + 0(h°°). 

4. Since X is compactly microlocalized, so are the operators Bi = X exp(— itA/h) and 
B 2 = Xexp{itA/h). Next, if Y 2 ,Y X G ^° mp , then 

Y 2 B 1 Y 1 = Y 2 X(exp(-itA/h)Y 1 exp(itA/h))exp(-itA/h); 

using our knowledge of the wavefront set of the operator in brackets, we see that this is 
0(h°°) if 

WF h (Y 2 ) n exp(tH ao ) WF h (Yi) = 0. 

Therefore, WFh(Bi) is contained in the graph of exp(tH ao ); similarly, WF^(B 2 ) is contained 
in the graph of exp(— tH ao ). Next, 

B X B 2 = X(exp(-itA/h)X exp{itA/h)); 

however, the operator in brackets is the identity microlocally near K, as X is the identity 
micro locally near K and K is invariant under exp(tH ao ). Therefore, BiB 2 is the identity 
microlocally near K. The intertwining property is proved in a similar fashion. 

5. We have 

d t (exp(tVX) exp(-tV)*) = exp(tVX)V(X - I) exp(-tV)* = 0(h°°) 
microlocally near K x K. □ 

Finally, we consider the special case a = 0; in other words, we study exp(itA), where A G 
^comp gj nce assoc i a ted canonical transformation is the identity, it is not unexpected 
that QxpiitA) is a pseudodifferential operator: 

Proposition 2.7. Let ai be the principal symbol of A. Then: 

1. expiitA) — I G \5/ c ° mp and the principal symbol of exp(itA) equals e ltai . Moreover, 
if A 1 = A 2 microlocally in some open set, then exp(itAi) = exp(itA 2 ) microlocally in the 
same set. 
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2. For any P G ^° mp ; we have the following asymptotic sum: 

exp(iL4)Pexp(-iL4) ~ 2^ ^ ^ — , 

i>o J ' 

where ad^ Q = [A, Q] for every Q. 

3. If U C T*M is connected and exp(iA) = I microlocally in U , then A = 2irl microlo- 
cally in U , where I is an integer constant. 

Proof. 1. We can find a family of pseudo-differential operators B t solving 

d t B t = iAB t + 0(h°°), B = I, 

by subsequently finding each member of the asymptotic decomposition of the full symbol 
ofB t . Then 

dt(exp(-itA)B t ) = 0(h°°); 
therefore, exp(itA) = B t + 0(h°°). The properties of B t can be verified directly. 

2. Follows directly from the equation 

d t (ex-p(itA)P exp(-itA)) = i &d A (exp(itA)P ex-p(-itA)). 

3. By calculating the principal symbol of exp(iA), we see that a\ has to be equal to 
2-71"/ in U for some constant / G Z. Subtracting this constant, we reduce to the case when 
A = 0(h). However, if A = 0(h N ) for some N > 1, then exp(iA) = I + iA + 0(h N+1 ); by 
induction, we get A = 0(h N ) microlocally in U for all N. □ 

2.4. Integrable systems. Assume that M is a two-dimensional manifold and Pi,p2 are 
two real- valued functions defined on an open set U dT*M such that: 

• {PiiPi} = 0; 

• for p = (pi,P2) '■ U — > M. 2 and each p G p(U), the set p _1 (p) is compact and 
connected. 

We call such p an integrable system. Note that if V C 1R 2 is open and intersects p(U), 
and F : V — > M 2 is a diffeomorphism onto its image, then F(p) is an integrable system on 

We say that an integrable system p : U — > 1R 2 is nondegenerate on U, if the differentials 
of pi and P2 are linearly independent everywhere on U. The following two propositions 
describe the normal form for nondegenerate integrable systems: 

Proposition 2.8. Assume that the integrable system p is nondegenerate on U . Then: 

1. For each p G p(U), the set p _1 (p) C T*M is a Lagrangian torus. Moreover, the 
family of diffeomorphisms 

<f>t = exp(tiP Pl + t 2 H P2 ), t = (*i,t 2 ) G M 2 , 
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defines a transitive action ofM? on p l (p). The kernel of this action is a rank two lattice 
depending smoothly on p; we call it the periodicity lattice (at p). 

2. For each p £ p{U), there exists a neighborhood V(po) an d a diffeomorphism F : 
V -»■ R 2 onto its image such that the nondegenerate integrable system F(p) has periodicity 
lattice 27rZ 2 at every point. Moreover, if the Hamiltonian flow of pi is periodic with minimal 
period 2ir, we can take the second component of F(p) to be p 2 . 

3. Assume that V C p(U) is open and connected and Fx,F% : V — > M 2 are two maps 
satisfying the conditions of part 2. Then there exist A £ GL(2,Z) and b £ IR 2 such that 
F 2 = A-F 1 + b. 

Proof. This is a version of Arnold-Liouville theorem; see |Du| Section 1] for the proof. □ 

Proposition 2.9. Assume that p : U — >■ M. 2 , p' : U' — > M. 2 , are nondegenerate integrable 
systems with periodicity lattices 2nZ 2 at every point; here U C T*M, U' C T*M' . Take 
p £ p(U) Dp' ([/')■ Then: 

1. There exists a symplectomorphism $ from a neighborhood o/p _1 (po) ^ n T*M onto a 
neighborhood of (p') _1 (Po) ^ n T*M' such that p = p' o $. 



2. $ is exact, as defined in Section 2.3, if and only if 



v S = I °' S i J = 1,2, 

1i 



where jj and^ are some fixed (2ir -periodic) Hamiltonian trajectories of pj on p 1 (po) and 
p'j on (p / ) _1 (p ); respectively. 

Proof. Part 1 again follows from Arnold-Liouville theorem. For part 2, we use that the 
closed 1-form a s — Q*a' s on a tubular neighborhood of p _1 (po) is exact if and only if its 
integral over each jj is zero. Since jj lie in p _1 (po) and the restriction of do~ s = u s to 
p _1 (po) is zero, we may shift jj to make both of them start at a fixed point (xo,£o) £ 
p _1 (p ). Similarly, we may assume that both start at $(xo,Co)- But in this case 7^ = 
$ o jj and 



[ a s - $V /S = / a S - [ 



which finishes the proof. □ 

Next, we establish normal form for one-dimensional Hamiltonian systems with one de- 
generate point. For that, consider IR 2 ,^ with the standard symplectic form d£ A dx, and 
define ( = (x 2 + £ 2 )/2; then ( has unique critical point at zero and its Hamiltonian flow is 
27r-periodic. 

Proposition 2.10. Assume thatp(x,^) is a real-valued function defined on an open subset 
of M 2 and for some A £ M., 
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• the set Ka = {p < A} is compact; 

• p has exactly one critical point (xo>£o) Ka, p(^o,£o) < A, and the Hessian of p 
at (x ,£o) 2S positive definite. 

Then there exists a smooth function F on the segment \p(xo, £o), A], with F' > everywhere 
and F(p(xo, £o)) = 0, and a symplectomorphism ^ from Ka onto the disc {( < F(A)} C IR 2 
such that F(p) = ( o Moreover, F'(p(x ,^ )) = (det V 2 p(x , £o))~ 1//2 - If P depends 
smoothly on some parameter Z , then F and \I/ can be chosen locally to depend smoothly on 
this parameter as well. 



Proof. Without loss of generality, we may assume that p{x ,^ ) = 0. Recall that in one 
dimension, symplectomorphisms are diffeomorphisms that preserve both area and orien- 
tation. By Morse lemma, there exists an orientation preserving diffeomorphism from a 
neighborhood of (xq, £o) onto a neighborhood of the origin such that p = (o9. Using the 
gradient flow of p, we can extend to a diffeomorphism from Ka to the disc {( < A} such 
that p = £ o 0. Let J be the Jacobian of _1 ; then the integral of J inside the disc {( < a} 
is a smooth function of a. Therefore, there exists unique function F smooth on [0, A] such 
that F' > everywhere, F(0) = 0, and the integral of J inside the disc {F(() < a}, that 
is, the area of _1 ({F(£) < a}) = {F(p) < a} C Ka, is equal to 2na. 

Let be a diffeomorphism from Ka onto {( < F(A)} such that F(p) = (oQ (constructed 
as in the previous paragraph, taking F(p) in place of p) and let J be the Jacobian of _1 . 
We know that for < a < F(A), the integral of J — 1 over {( < a} is equal to 0. Introduce 
polar coordinates (r, ip); then there exists a smooth function ip such that J = 1 + d^ip (see 



Proposition 3.7). The transformation 

^ : (r, ip) h-> (r,<p + ip) 

is a diffeomorphism from {( < F(A)} to itself and has Jacobian J; it remains to put 
$ = $o0. To compute F'(p (x ,^ )), we can compare the Hessians of F(p) and £ o $ at 
(zo,6o)- 

The function F is uniquely determined by p and thus will depend smoothly on Z. As 
for \]/, we first note that was constructed using Morse lemma and thus can be chosen 
locally to depend smoothly on Z (see for example |EvZw[ Proof of Theorem 3.15]). Next, 
we can fix ip by requ iring that it integrates to zero over each circle centered at the origin 



see Proposition 3.7); then ip, and thus will depend smoothly on Z. □ 



3. Angular problem 
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3.1. Outline of the proof. Consider the semiclassical differential operators (using the 
notation of ( 1.11 )) 



Px(u,D;h) = h 2 P e (co 
(1 + a) 2 



sin 9 



(hD e )(A e sin9 ■ hD e 



A sin^ 9 



, n y a(uj + ihv) sin 2 9 — hD^) 2 



P 2 (h) = hD^ 



on the sphere § 2 . Then (co, A, k) is a pole of Rg if and only if (A + ihfi, k) lies in the joint 
spectrum of the operators (Pi,P 2 ) (see Definition A.l). For a = 0, Pi is the Laplace- 
Beltrami operator on the round sphere (multiplied by —h 2 ); therefore, the joint spectrum 
of (Pi,P 2 ) is given by the spherical harmonics (1(1 + h),k), k,l G hZ, \k\ < I (see for 
example |Tay[ Section 8.4]). In the end of this subsection, we give a short description of 
which parts of the angular problem are simplified for a = 0. For general small a, we will 
prove that the joint spectrum is characterized by the following 



Proposition 3.1. Let u,u satisfy (1.14); we suppress dependence of the operators and 
symbols on these parameters. Consider 

K = {(A, k) | C e l <\<C e , A > (1 + a) 2 {k - au) 2 } C IR 2 , 

K ± = {(A, k) G K | (1 + a)(k - auj) = ±Vl}. 

Then there exist functions G±(X, k; h) such that: 

1. G± is a complex valued classical symbol in h, smooth in a fixed neighborhood of K . 
For (A, k) near K and \Jx\ < Cg, we can define G±(X + ihp,, k) by means of an asymptotic 
(analytic) Taylor series for G± at (A, k). 



2. For a = 0, G±(\, k; h) = -h/2 + yj A + h 2 /A T k. 

3. G„{\,k;h)-G+{\,k;h) =2~k. 

4- Let F± be the principal symbol of G±. Then F± is real-valued, d^F± > and ^d~ k F± > 
on K , and F±\j£ = 0. 

5. For h small enough, the set of elements (A + ihjl, k) of the joint spectrum of (Pi, P?) 
satisfying (1.14) lies within 0(h) of K and coincides modulo 0(h°°) with the set of solutions 
to the quantization conditions 

k E hZ, G ± (X + ihfl,k) e KN; 

here N is the set of nonnegative integers. Note that the conditions G + G Z and G_ G Z 
are equivalent; however, we also require that both G + and G_ be nonnegative. Moreover, 
the corresponding joint eigenspaces are one-dimensional. 
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Figure 2. The joint spectrum and the set K. 



1.6 follows from the proof of Proposition 3.1 In fact, the symbol T (I, u, u, k; h) 



Proposition 

is defined as the solution A + ihjl to the equation 



G+(A + ihfi, k, uj + ihv\ h) = I — k\ 



this proves part (1) of Definitio n |1.1 The resolvent estimates are an immediate corollary 
of the ones stated in Proposition 3.8 below. The decomposition of T§ at a = follows from 



Proposition 3.4 



We now give the schema of the proof of Proposition 3.1 Let pjo be the principal symbol 
of Pj] note that both p w and P20 are real-valued; also, define p = (j»io,P2o) : T*S 2 — > M. 2 . 
In Section 3.2, we construct the principal parts F± of the quantization symbols globally in 
K, and show that the intersection of the image of p with {C^ 1 < A < Cg} is exactly K. 



Using the theory of integrable systems described in Section 2.4, we then construct local 
symplectomorphisms conjugating (F±(p) 
where Ai = 
(x,y), and 



P20) away from K T to the system ((, rj) on T*A4, 
is called the model space, (£, 77) are the momenta corresponding to 



„2 + {:2 



Note that the integrable system ((, rj) is nondegenerate on {( > 0} with periodicity lattice 
2ttZ 2 , and d( = on {C = 0}. 

Next, we take (Ao, k^) e K and show that joint eigenvalues in a certain /^-independent 
neighborhood of this point are given by a quantization condition. For this, we first use 



Egorov's theorem and the symplectomorphisms constructed in Section 3^ to conjugate 
P\,P2 microlocally near p _1 (Ao,fco) to some pseudodifferential operators Qi,Q2 on M.. 
The principal symbols of Qj are real-valued functions of (C,^) only; in Section 3.3, we use 
Moser averaging to further conjugate Qi, Q 2 by elliptic pseudodifferential operators so that 
the full symbols of Qj depend only on (C,^)- In Section 3.4, we use spectral theory to 



construct a local Grushin problem for (Qi,Q2), which we can conjugate back to a local 
Grushin problem for (Pi, P2); then, we can apply the results of Appendix |A] to obtain local 



ASYMPTOTIC DISTRIBUTION OF QNMS FOR KERR-DE SITTER 



20 



quantization conditions (Proposition 3.8). To pass from these local conditions to the global 
one, we use 

Proposition 3.2. Assume that Gj(X,k;h) are two complex-valued classical symbols in h 
defined in some open set U C M 2 , their principal symbols are both equal to some real-valued 
F(X,k), with d~ x F ^ everywhere and {F > 0} convex, and solution sets to quantization 
conditions 

k G hZ, Gj{~X + ihfi,k) G hN 

in the region (A, k) G U, jl = 0(1) coincide modulo 0(h°°). Then G\ — G 2 = hi + 0(h°°) 
on {F > 0} for some constant I G Z. Moreover, if {F = 0} D U 7^ 0, then 1 = 0. 

Proof. Assume that (X±,ki) G {F > 0}. Then for every h, there is a solution (X(h) + 
ihfi(h),k(h)) to the quantization conditions within 0(h) of (X±,ki); we know that 

Gj(X(h) + ihjl(h), k(h)) ehZ + 0(h°°), j = 1, 2, 

and thus (G x - G 2 )(X(h) + ihp,(h),k(h)) = hl(h) + 0(h°°), for some 1(h) G Z. Since 
G\ — G2 = O(h) in particular in C 1 , we have 

\(G 1 -G 2 )(X 1 ,k 1 )-hl(h)\=0(h 2 ). 

Therefore, 1(h) is constant for h small enough and it is equal to the difference of subprincipal 
symbols of Gi and G 2 at (Xi,kx). It follows that 1(h) is independent of (Ai,&i); we can 
subtract it from one of the symbols to reduce to the case when G\ — G 2 = 0(h 2 ). The 
analysis in the beginning of this proof then shows that 

\\G\ — G2||c(f>o) = 0(h\\Gi — G 2 ||c 1 (f>o) + h°°). 

Arguing by induction, we get G± — G 2 = 0(h N ) for all N. The last statement follows 
directly by taking solutions to the quantization conditions with Gj = and requiring that 
they satisfy the quantization conditions G%-j > 0. □ 

We can now cover K by a finite family of op en sets, on each of which there exists a local 



quantization condition. Using Proposition 3.2 and starting from K±, we can modify the lo- 
cal quantization conditions and piece them together to get unique (modulo h°°) global G±. 
The joint spectrum of (Pi, P 2 ) in a neighborhood of K is then given by the global quanti- 



zation condition; the joint spectrum outside of this neighborhood, but satisfying (1.14), is 



empty by part 2 of Proposition 3.8 



Also, the principal symbol of G_ — G + is equal to 2k; therefore, G_ — G + — 2k is equal to 
Ih for some fixed I G Z. However, G± depend smoothly on a and thus it is enough to prove 
that I = for a = 0; in the latter case, the symbols G± are computed explicitly from the 
spectrum of Laplacian on the round sphere. (Without such a reference point, one would 
need to analyse the subprincipal symbols of G± using the Maslov index.) This finishes the 



proof of Proposition 3.1 
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Finally, let us outline the argument in the special case a = and indicate which parts of 
the construction are simplified. The formulas below are not used in the general argument; 
we provide them for the reader's convenience. The principal symbol pio of Pi is just the 
square of the norm on T*§ 2 generated by the round metric: 



£ 2 

PlO = 60 + 



sin 2 6 



The set p 1 (A, k) consists of all cotangent vectors with length v A and momentum k; 
therefore 

(1) for A < k 2 (corresponding to the complement of K), the set p~ 1 (X,k) is empty; 

(2) for k = ±\/a (corresponding to K±), the set p _1 (A,fc) is a circle, consisting of 

covectors tangent to the equator with length vA and direction determined by the 
choice of sign; 

(3) for A > k 2 (corresponding to the interior of K), the set p _1 (A, k) is a Liouville torus. 



can 



3.5 



The principal parts F± of the quantization symbols, constructed in Proposition 3.4 
be computed explicitly: F± = VT =F k (see the proof of part 2 of this Proposition). T hen 
F± {Ci 7 )) = (C =t v) 2 - For > 0) the canonical transformation $-1- from Proposition 
can be taken in the form 

= ((2pio) 1/2 ( v / ^ ± ^Y 1 ' 2 c os9, ip + G, -2 1 / 2 (VP^ ± U~ 1/2 sin^, ' 



here (x, y, rj) are coordinates on T*A4, with Ai = R x x §^ the model space. The function 
G : T*§ 2 -> S 1 here is given by 

pT Q ± £ v ) cos G = pK? sin 6> ± 4^, {\/p~io ± C v ) sin G = =f cos 6£ e . 

sm (7 



In fact, the maps $-t defined in (3.1 ) extend smoothly to the poles {sin 6 = 0} of the sphere 



and satisfy the conditions of Proposition 3.5 on the complement of the opposite equator 

{6 = tt/2, & = 0, ^ = TV 7 !}. 

One can then conjugate the operators Pi, P 2 to some model operators Qi, Q 2 as in Propo- 
3.6 To bring the subprincipal terms in Qj to normal form, one still needs Moser 



sit ion 



averaging. Once the normal form of Proposition |3.6 is obtai ned, i t is possible to use the 



ellipticity of p — (A, k) away from p l (X, k) (as in Proposition A. 4) and spectral theory to 



obtain the quantization condition. The Grushin problem construction of Section |3.4| and 



Appendix |A.1| is not needed, as the operator Pi is self-adjoint. 



3.2. Hamiltonian flow. Let (9, if) be the spherical coordinates on § 2 and let (£0,^) be 
the corresponding momenta. Note that £0 is defined away from the poles {sin# = 0}, 
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while £ v is well-defined and smooth on the whole T*§ 2 . In the (9, ip, £g, ^) coordinates, the 
principal symbols of P 2 and Pi are p 2 o = £ v and 



l + af 
A,sin 2 VSv 



Since pio does not depend on <p, we have 

{PlO,P2o} = 0. 



We would like to apply the results of Section 2.4 on integrable Hamiltonian systems to 
establish a normal form for p = (pio,P2o)- First of all, we study the points where the 
integrable system p is degenerate: 

Proposition 3.3. For a small enough, 

1. For Cq 1 < A < Cg, the set p _1 (A, k) is nonempty if and only if (A, k) E K . 

2. The integrable system p is nondegenerate on p _1 (if) ; except at the equators 

E±(X) = {9 = vr/2, ie = 0, (1 + a){i v - au) = iv 7 !} C T*§ 2 , C e x <~\<C Q . 
Moreover, p w = A on E±(A) and the union of all E±(\) is equal to p -1 (fC|-). Also, 



dp 10 = ±2(1 + a)V\dp 20 on E±(\). (3.2) 



Proof. We can verify directly the statements above for a = 0, and also (3.2) for all a. 
Then part 2 follows for small a by a perturbation argument; part 1 follows from part 2 by 
studying the extremum problem for £ v restricted to {p w = A}. □ 

Next, we construct the principal parts F± of the quantization symbols globally: 

Proposition 3.4. For a small enough, 

1. There exist unique smooth real-valued functions F±(\,k) on K such that F±\f£ =0 

and (F±(p),p 2 o) is a nondegenerate completely integrable system on p^ x (K \ (K + U K-)) 
with periodicity lattice 2ttZ* 2 . 

2- d^F± > 0, Td%F± > 0, and F-(X,k) — F + (X,k) = 2k on K. In particular, one 
can define the inverse P± 1 (C)^) of F± in the A variable, with k as a parameter. Also, 
F± = V\Tkfora = 0and d~ x F ± = ±(2k)~ 1 + 0(a 2 ) on K±. 

3. If (A, k) G K \ (K + U KJ) and j± are some (2tt -periodic) trajectories of F±(p) on 
p _1 (A, k), then 

a s = 27iF±(\,k). (3.3) 

7± 
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Proof. 1. We first construct F + in a neighborhood of K + . In fact, we take small > 
and define F + on the set 

K £ = {(A, jfe) | jfe > e k , (1 + a) 2 (A; - auf < A < C e }. 

We will pick small enough so that K + C iC e ; note, however, that does not lie in X. 
Moreover, we will construct a symplectomorphism $ from p _1 (il" e ) onto a subset of T*Ai 
such that -F+(p) ° = C an d P20 ^ 1 — 7 7- 

Note that ^ = on the poles of the sphere {sin 9 = 0}; therefore, (0, y?, £0, £^) is a 
symplectic system of coordinates near p _1 (il" e ). Next, fix £^ > e fc and consider p 10 as 
a function of (0,£e)j then for a small enough, this function has a unique critical point 
(0,0) on the compact set {pio(-, — the Hessian at this point is positive definite. 
Indeed, it is enough to verify these statements for a = and check that dePw = (%,Pio = 
for (0, £0) = (7r/2,0) and small a. Now, we may apply Proposition 2.10 to the function 
{Piq("> ') £<p)} an d obtain a function F + (X] k) on K £ such that = and d^F + > and 

a mapping 

* : (0, H- (^(0,&,^),^(0,&,^)) 

that defines a family of symplectomorphisms (0, £0) 1— )• (\I/ X , \]>g), depending smoothly on 
the parameter and 

Now, define $ : (0,^,6?,^) ($ x , $^) G T*.M by 

Here G(x,£;£ v ) is some smooth function. For $ to be a symplectomorphism, G should 
satisfy 

Since vary in a disc, this system has a solution if and only if 

= d^ x } + ¥,} = VI/,.}; 

this is true since {^g, *& x } — 1- The defining properties of F + now follow from the cor- 
responding properties of the integrable system uniqueness follows from part 3 of 
Proposition 



2.8 



and the condition F + \g = 0. 



Now, by part 2 of Proposition 2.8 and Proposition 3.3, for each (Ao, fco) £ ^\ (-K+UJf_), 
there exists a smooth function -F(A, k) defined in a neighborhood of (Ao, ho) such that 
d\F 7^ and (F(p),p 2 o) has periodicity lattice 27rZ 2 ; moreover, part 3 of Proposition 2.8 



describes all possible F. Then we can cover K \ K £ by a finite set of the neighborhoods 
above and modify the resulting functions F and piece them together, to uniquely extend 
the function F + constructed above from K e to K \ . (Here we use that K \ K_ is simply 
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connected.) Similarly, we construct F_ on K\K + . (The fact that F± is smooth at K T will 
follow from smoothness of F T at K T and the identity F_ — F + = 2k.) 

2. We can verify the formulas for F± for a = explicitly, using the fact that the 

Hamiltonian flow of y/% is 27r-periodic in this case. The first two identities now follow 
immediately. As for the third one, we know by part 3 of Proposition 2JS and the case a = 
that F_ — F + = 2k + c for some constant c; we can then show that c = using part 3 of 
this proposition. Finally, d^F± \ g can be computed using Proposition 



2.10 



3. First, assume that (A, k) G K £ and let $ be the symplectomorphism constructed in 
part 1. Then 

$ o 7 + = {( = F + (X, k), r) — k, ip = const} 
is a circle. Let D + be the preimage under $ of the disc with boundary $ o 7 + ; then 



a 



7+ 



to 



M 



2nF + (\,k). 



We see that (3.3) holds for F + near K + ; similarly, it holds for F_ near K_. It now suffices 
to show that for each (Ao, ko) G K \ (K + UK J), there exists a neighborhood V(Xo, k ) such 
that if (Xj,kj) G V, j = 1,2, and 7^ are some (27r-periodic) Hamiltonian trajectories of 
F±(p) on p _1 (Aj,fcj), then 



a s = 2vr(F ± (A 2 ,fc 2 )-F ± (A 1 ,fc 1 )). 



(3.4) 



Ti 



In particular, if (3.3) holds for one point of V, it holds on the whole V. One way to 



prove (3.4) is to use part 1 of Proposition 2.9 to conjugate (F±(p),p 2 o) to the system 
(£x)£jj) on the torus T x<y and note that the left-hand side of ( |3.4 ) is the integral of the 
symplectic form over a certain submanifold bounded by 71,72; therefore, it is the same for 
the conjugated system, where it can be computed explicitly. □ 

Finally, we construct local symplectomorphisms conjugating (F±(p),p 2 o) to (C,v) : 

Proposition 3.5. For each (Ao, ko) & K \ K T , there exists an exact symplectomorphism 
$-t from a neighborhood of p~ 1 (Ao, ko) in T*S 2 onto a neighborhood of 

J^M = {C = F±(\ , ko), 7] = ko] 

in T*A4 such that 



Pio o = F^iCr]), p 2 o o 



2.9 



applied 



Proof. The existence of $± away from K± follows from part 1 of Proposition 
to the systems (F±(p),p 2 o) and (Cj^)j near K±, these symplectomorphisms have been 
constructed in the proof of part 1 of Proposition 3.4 Exactness follows by part 2 of 



Proposition 2.9 (which still applies in the degenerate case); the equality of classical actions 
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over the flows of i*±(p) and ( follows from part 3 of Proposition 3.4, while the classical 



actions over the flows of both p 2 o on p 1 (Ao, ko) and 77 on Am are both equal to 2nko- □ 

3.3. Moser averaging. Fix (Ao, ko) G K \ K T , take small e > 0, and define (suppressing 
the dependence on the choice of the sign) 

A = p-^Ao, fco), Co = F±(A , h), A° M = {( = Co, V = ~k }c T*M, 

V £ = {(A, k) I |F±(A, k) - Co| < e, |fc - fco| < £> C K 2 , 

V% = {IC - Col < e, \n- fcol < e} C T*M; 

then and are compact neighborhoods of (Aq, ko) and A^, respectively. Here the 



functions F± are as in Proposition 3.4 Let $-1- be the symplectomorphism constructed in 
Proposition 3.5 we know that for e small enough, 3>±(p _1 (V e )) = Vj^. In this subsection, 



we prove 

Proposition 3.6. For (Aq, ko) & K \ K± and e > small enough, there exists a pair of 



operators (Bi,B 2 ) quantizing $± near p 1 (V e ) x in the sense of Proposition 2.5 and 
operators Qi,Q 2 G \I/™ mp (.M) such that: 

1. Pi and P 2 are intertwined with Qi and Q 2 , respectively, via (B 1 ,B 2 ), nearp^iV 6 ) x 
Vjw- P follows immediately that the principal symbols of Q± and Q 2 are i 7 ^ 1 77) and rj, 
respectively, nearVj^. 

2. Q 2 = hD y and the full symbol of Qi is a function of (Crj), microlocally near Vj^. 
Here we use Weyl quantization on M., inherited from the covering space M. 2 . 



First of all, we use Proposition 2.5 to find some (Bx,B 2 ) quantizing $ ± and Qi,Q 2 
intertwined with Pi and P 2 by (Bi,B 2 ). Then we will find a couple of operators 1,7 6 
*7 P (M) such that Y = X" 1 near V e M and the operators Q[ = XQ X Y,Q' 2 = XQ 2 Y 



satisfy part 2 of Proposition |3.6 This is the content of this subsection and will be done in 



several conjugations by pseudodifferential operators using Moser averaging technique. We 



can then change Bi,B 2 following the remark after Proposition |2.5| so that Pi and P 2 are 
intertwined with Q[ and Q 2 , which finishes the proof. 

The averaging construction is based on the following 

Proposition 3.7. Assume that the functions po, /o, g G C°°(V^) are given by one of the 
following: 

(1) po = /o = T] and g is arbitrary; 

(2) Po = ( and f = fo((,v) ^ s smooth in Vf^, with d^f 7^ everywhere, and g is 
independent ofy. 

Define 

1 f 2w 

(9) = 7T / 9 e w(tH po ) dt. 
^ Jo 
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Then there exists unique b £ C°°(V^) such that (b) = and 

9= (g) + {fo,b}. 

Moreover, in case (2) b is independent of y. 

Proof. We only consider case (2); case (1) is proven in a similar fashion. First of all, if b is 
y-independent, then {/ , b} = d^f -{(, b} = {(, d^f -b}; therefore, without loss of generality 
we may assume that fo — C- The existence and uniqueness of b now follows immediately if 
we treat y, r\ as parameters and consider polar coordinates in the (x,£) variables. To show 
that b is smooth at £ = (in case (o < e) , let z = x + i^ and decompose g — (g) into an 
asymptotic sum of the terms z^z k with j, k > 0, j ^ k, and coefficients smooth in (y,rj); 
the term in b corresponding to z^z k is z^z k j{i(k — j)). □ 

Henceforth in this subsection we will work with the operators Qj on the level of their full 
symbols, microlocally in a neighborhood of V%^- (The operators X and Y will then be given 
by the product of all operators used in conjugations below, multiplied by an appropriate 
cutoff.) Denote by qj the full symbol of Qj. We argue in three steps, following in part jHiSj 
Section 3]. 

Step 1: Use Moser averaging to make q2 independent of y. 

Assume that q 2 is independent of y modulo 0(h n+1 ) for some n > 0; more precisely, 

n 

q 2 = J2 hj( l*Ax,t;,v) + h n+1 r n (x,y,C,v) + 0(h n+2 ). 

Take some B e \I/J.° C with principal symbol b and consider the conjugated operator 

Q' 2 = exp(ih n B)Q 2 exp(-ih n B). 



Here exp(±ih n B) £ W°\ c are well-defined by Proposition 2.7 and inverse to each other; 



using the same proposition, we see that the full symbol of Q' 2 is 

n 

h j q2, 3 {^ & V) + h n+1 (r n - b}) + 0(h n+2 ). 

3=0 



If we choose b as in Proposition 3.7 1), then r n — {77, b} = (r n ) is a function of (x, rj) only; 
thus, the full symbol of Q' 2 is independent of y modulo 0{h n+2 ). Arguing by induction and 
taking the asymptotic product of the resulting sequence of exponentials, we make the full 
symbol of Q 2 independent of y. 

Step 2: Use our knowledge of the spectrum of P 2 to make q 2 = 7]. 
First of all, we claim that 

exp(27riQ 2 //i) = I (3.5) 
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microlocally near 

V M x V M- For that > we wil1 use Proposition Y2M Let X G ^ omp (§ 2 ) have 
real- valued principal symbol, be microlocalized in a small neighborhood of p _1 (V^ £ ), but 
equal to the identity microlocally near this set. Consider 

V t = exp(ttQ 2 /h)B 1 exp(-it(P 2 X)/h)B 2 . 

We see that 

hDtPt = exp(itQ 2 /h)(Q 2 B 1 - B X P 2 X) exp(-it(P 2 X)/h)B 2 
vanishes microlocally near Vf^ x V^; integrating between and 2-7T and using part 5 of 



we get (3.5 ). 



Proposition 2.6 to show that exp(— 2iri(P 2 X) / h) = I microlocally near p 1 (V e ) x p 1 (V e ), 



Now, let Xm G \I/™ mp (.M) be equal to the identity microlocally near WFh(Q 2 )', since 
the full symbol of Q 2 is independent of y, we have [Q 2 , {hD y )X M } = 0(h°°). Therefore, by 
parts 2 and 5 of Proposition |2.6| 

exp(27ri(Q 2 — (hD y )XM)/h) = exp(— 2niD y XM) exp(2niQ 2 /h) = I 

microlocally near x Vj^. However, R = h~ 1 (Q 2 — (hD y )XM) £ near and 



thus the left-hand side exp(2iriR) is pseudodifferential; by part 3 of Proposition 2.7 , we get 
R = I for some constant I G Z and therefore 

Q 2 = hD y + hi 

microlocally near Vf^. It remains to conjugate Q 2 by e tly to get q 2 = rj. 

Step 3: Use Moser averaging again to make q\ a function of ((, 77), while preserving q 2 = rj. 

Recall that [Pi,P 2 ] = 0; therefore, [Qi,Q 2 ] = (microlocally near V^). Since q 2 = r], 
this means that qi is independent of y. We now repeat the argument of Step 1, using 



Proposition 3.7^2) with f = i ? ± 1 (C ) ? 7)- The function b at each step is independent of y\ 
thus, we can take [B, hD y ] = 0(h°°). But in that case, conjugation by exp(ih n B) does not 
change Q 2 ; the symbol of the conjugated Qi is still independent of y. Finally, (r n ) is a 
function of (C,^); therefore, q\ after conjugation will also be a function of (C, 7 ?)- 

3.4. Construction of the Grushin problem. In this subsection, we establish a local 
quantization condition: 

Proposition 3.8. 1. Assume that (Ao, fco) £ K\K T and V £ is the neighborhood o/(Ao,fco) 
introduced in the beginning of Section 3.3[ Then for e > small enough, there exists a 



classical symbol G±(\, k; h) on V £ with principal symbol F± and such that for l:6Z, the 
poles A + ihpL of Rg(u, A, k) with (X,k) G V £ and < C g are simple with polynomial 
resolvent estimate L 2 — > L 2 , in the sense of Definition \ l and coincide modulo 0(h°°) 
with the solution set of the quantization condition 

G±(\ + ihp,,k;h) G hN. (3.6) 
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2. Assume that (X ,k ) satisfies (1.14), but does not lie in K. Then there exists a 
neighborhood V{\q, Hq) such that for h small enough, there are no elements (A + ihfi, k) 
of the joint spectrum of Pi,P 2 with (X,k) G V and \fi\ < Ce, and Rg(u,X,k) is bounded 
L 2 ^L 2 byO{h 2 ). 

To prove part 1, we will use the microlocal conjugation c onst ructed above. Let (Aq, ko) G 



K \ K T and e > 0, B\, B 2 , Qi, Q 2 be given by Proposition 3.6 Consider the operators 



T l = l -((hD x ) 2 + x 2 )-^ T 2 = hD y 

on their full symbols are ( — h/2 and rj, respectively. We know that Ti and T 2 commute; 
the joint spectrum of T 1; T 2 is h(N x Z). Therefore, for any bounded function / on IR 2 , we 
can define f(T 1 ,T 2 ) by means of spectral theory; this is a bounded operator on L 2 (Ai). 

Proposition 3.9. 1. For f G C^°(IR 2 ) ; the operator f(Ti,T 2 ) is pseudodifferential; more- 
over, f(Tt,T 2 ) G ^ mp (M) and WF h (/(Ti,r 2 )) C {(C,v) e supp/}. The full symbol of 
f(Ti,T 2 ) in the Weyl quantization is a function of ( and 77 only; the principal symbol is 

f(C,v)- 

2. Assume that d G hN, rji G KL. Let u be the L 2 normalized joint eigenf unction of 
(Ti,T 2 ) with eigenvalue (Ci ? ^7i) ■ Then u is compactly microlocalized and 

WF fc (w) C {C = Ci, V = Vi}- 



3. Assume that the function f((,r];h) is Borel measurable, has support contained in a 
compact h-independent subset Kf ofM 2 , and 

max{\f((,ri; h)\ \ ( G hN, 7] G KL} < ChT r 

for some r > 0. Then the operator f(Ti,T 2 ;h) is compactly microlocalized, its wavefront 
set is contained in the square of {(C,v) £ Kf}, and the operator norm of f(Ti,T 2 ; h) is 
0{h~ r ). 

Proof. For part 1, we can show that the operator f(Ti,T 2 ) is pseudodifferential by means 
of Helffer-Sjostrand formula in calculus of several commuting pseudodifferential operators; 
see for example |DiSj , Chapter 8]. This also gives information on the principal symbol 



and the wavefront set of this operator. To show that the full symbol of f(Ti,T 2 ) depends 
only on ((,r)), note that if A G ^ loc (Ai) and a is its full symbol in the Weyl quantization, 
then the full symbol of [A, Ti] in the Weyl quantization is — ih{a,(}] similarly, the full 
symbol of [A, T 2 ] in the Weyl quantization is —ih{a, rf} (see for example |Sj2[ discussion 



before (1.11)]). Since [/(Ti, T 2 ), Tj] = 0, the full symbol of f(Ti,T 2 ) Poisson commutes 
with C and rj. 

To show part 2, we take x(C> r ?) C^°(1R 2 ) equal to 1 near (£i,?7i); then u = x(Ti,T 2 )u. 
Similarly, to show part 3, we take x equal to 1 near Kf] then the L 2 operator norm of 
/(Ti, T 2 ) can be estimated easily and f(T 1 ,T 2 )= x(T 1; T 2 )f(T 1 , T 2 )x(Ti, T 2 ). □ 
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Now, recall that by Proposition 3.6, the full symbol of Q% in the Weyl quantization is a 
function of ((, r/) near V^; therefore, we can find a compactly supported symbol G±(X, k\ h) 
such that the principal symbol of G± near V e is F± and 

Qi = G ± 1 (T ll T 2 ;h) 

microlocally near Vj^, where G± (£, k; h) is the inverse of G± in the A variable. Recall also 
that Q 2 = T 2 microlocally near Vj^. Multiplying Qi,Q 2 by an appropriate cutoff, which 
is a function of T 1; T 2 , we can assume that Q±,Q 2 are functions of T 1; T 2 modulo h 00 ^/^ 00 . 
We can now construct a local Grushin problem for Qi,Q 2 : 

Proposition 3.10. Let (Ai,^i) G V s and < Cq. 



1. Assume that (Ai + ihjji,ki) satisfies (3.6), with = G±(Ai +ihfli,ki) G hN. Then 

are 



there exist operators Ai,A 2 ,Si,S 2 such that conditions (L1)-(L5) of Appendix A. 2 



satisfied, with r — 1, (Pi, P 2 ) replaced by (Qi — X\ — ihjli,Q 2 — ki), K — {£ = £i, r\ — k\\, 
and 

A 1 {Q 1 - Ax - i/i/ii) + ^ 2 (Q 2 -h)=I- (3.7) 
microlocally near x . 
U. Fix 5 > and assume that 

\{Xi + ihfiuk) -{Gl l {C,v),v)\ >$h, (ehN, v e hZ. 



Then there exist operators Ai,A 2 such that the conditions (L1)-(L2) of Appendix A. 2 are 
satisfied, with r = 1, (Pi, P 2 ) replaced by (Qi — Ai — ihfl 1 ,Q 2 — k\), K — {£ = £i> r) — k\}, 
and 

Ai(Qi - Ai - ih[i x ) + A 2 (Q 2 -h)=I 
microlocally near x . 

Proof. 1. Let S\ : C — > L 2 (M) and S2 : L 2 (A4) — >■ C be the inclusion and the orthogonal 
projection onto, respectively, the unit joint eigenfunction of (T 1; T 2 ) with eigenvalue (Ci, ki). 



The properties (L3) and (L4) now follow from part 2 of Proposition 3.9 



Next, we use a partition of unity on the circle to construct the functions X11X2 with the 
following properties: 

• Xj C* 00 ^ 2 \ 0) is positively homogeneous of degree 0; 

• Xj > an d Xi + X2 — 1 everywhere on IR 2 \ 0; 

• Xj(s u s 2 ) =0for \sj\ < |s 3 -j|/2. 

It follows that 

l^d^ilM^I^I^Cd^il + l^l)" 1 - (3-8) 
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Take xidv) £ Co° supported in a small neighborhood of Vj^, while equal to 1 near Vj 
define the functions f\, f 2 as follows: 



Mi 



x(C, v)Xi(\G± 1 (C, ri; h) - Ai - ih^l 2 , \y - h\ 

G^iCr]; h)-Xi- ihfit 
x(C,v)X2(\G^}(C,V, h ) - Ai -ihjx^ 2 , \rj - h\ 



for (£, 77) 7^ (d, fci); we put /j(Ci, &i) = 0. We now take Aj = fj(Ti, T 2 ; h). Noticing that 

\Gg\C, m h)-\ x \ + \n-ki\> h/c, ((, v) e &(N xZ)n su PP x \ (Ci, fei), 

and using Proposition 3.9 and (3.8), we get that Aj are compactly microlocalized and 
ll^ill = Moreover, if x((,rj) is equal to 1 near (d, ki), then (1 — x)fj are smooth 

symbols; then, A" 



[1 — x)(Ti,T2)Aj belongs to \j>^° mp \yy p ar t \ Q f Proposition 



3.9 



and 
nave 



Aj = x(Tx,T 2 )Aj is microlocalized in the Cartesian square of {((,?]) G suppx}; we 
established property (LI), with r — 1. The properties (L2), (L5), and (3.7) are easy to 
verify, given that all the operators of interest are functions of T 1; T 2 . 

2. This is proved similarly to part 1. □ 

Finally, we conjugate the operators of the previous proposition by Bi,B 2 to get a local 
Grushin problem for P 1; P 2 and obtain information about the joint spectrum: 



Proof of Propositon 3.8 1. Assume first that \\,k\,pL\ satisfy the conditions of part 1 of 
Proposition 3.10; let Ai, A 2 , Si, S 2 be the operators constructed there. Recall that A,- are 
microlocalized in a small neighborhood of V%^- Then the operators 

Aj = B 2 AjB>i, Si = B 2 Si, S 2 = S 2 Bi, 

together with Pi — Ai — ihfii, P 2 — k\ in place of Pi, P 2 satisfy the conditions (L1)-(L5) of 
Appendix A. 2 with K = p _1 (Ai, k\) and 



Ai{Pi - Ax - ihfii) + A 2 (P 2 -ki) = I- SiS 2 

microlocally near p _1 (V e ). Moreover, Pi — Ai — ihpLi, P 2 — ki satisfy conditions (E1)-(E2) 
of Appendix A. 2 and the set where both their principal symbols vanish is exactly K. We 
can now apply part 2 of Proposition A. 4 to show that for h small enough and some 6 > 0, 
independent of h, \i,k\, there is exactly one element of the joint spectrum of (Pi, P 2 ) in the 
ball of radius Sh centered at (Ai + ihjli, ki), and this point is within 0(h°°) of the center 
of the ball. 

Now, we assume that (A + ihfl, k) satisfies the conditions of part 2 of Proposition 



3.10 



with 8 specified in the previous paragraph. Then we can argue as above, using part 1 o 
Proposition |A.4 , to show that this point does not lie in the joint spectrum for h small 
enough. 
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Since every point (A + ihp,, k) such that (A, k) 6 V s and < Cq is covered by one of 
the two cases above, we have established that the angular poles in the indicated region 
coincide modulo 0(h°°) with the set of solutions to the quantization condition. Moreover, 



Proposition A.3 together with the construction of a global Grushin problem from a local 



one carried out in the proof of Proposition A. 4 provides the resolvent estimates required in 
Definition 11.11 



2. The set p (Ao, ko) is empty by Proposition 
T = (P x - A - z/i//)*(Pi - A 



3.3 therefore, the operator 
ihp) + (P 2 - k) 2 



2 ) for (A, k) close to (Ao, k ) and /i bounded; therefore, for h 
= 0(1). The absense of joint spectrum and resolvent estimate 

□ 



is elliptic in the class \I /2 ( 
small enough, ||T _1 || L 2^ L 2 

follow immediately if we notice that the restriction of T to V' k is h 4 (Pe — X)*(Pe — A) 



4. Radial problem 



4.1. Trapping. In |Dyl[ Section 4], we use a Regge- Wheeler change of variables r — > x, 
under which and after an appropriate rescaling the radial operator becomes (using the 



notation of (1.11)) 



Px(h) = h 2 D 2 + V(x, u>, v, A, fi, k; h), 
V(x; h) = (A + ihp)A r - (1 + a) 2 ((r 2 + a 2 )(uj + ihv) - akf 



(note the difference in notation with |Dyl[ Section 7]). Let V(x; h) = Vq(x) + hV\(x) + 
h 2 V2(x), where 

V (x) = AA r - (1 + a) 2 ((r 2 + a 2 )u - ak) 2 

is the semiclassical principal part of V(x); note that Vq is real- valued and for 1 < £j < 2 
and a small enough, Vo(±oo) < 0. Now, |Dyl Proposition 7.4] establishes an arbitrarily 
large strip free of radial poles in the nontrapping cases; therefore, the only radial poles in 



the region (1.12) appear in case (3) of |Dyl Proposition 7.3]. Using the proof of the latter 
proposition, we may assume that: 

• | A — Ao(o), ^)| < £ r , where Ag 1 is the value of the function 

A 

F v (r;u,k) 



(l + a) 2 ((r 2 + a 2 )u-aky 



at its only maximum point. Under the assumptions ( 1.12[ ), l/C < Ao < C for some 
constant C; 

V , as a function of x, has unique global maximum x , |Vo(x )| < £^ and Vq(x) < 

for \x — x \ < e r ; 

Vq(x) < —el for \x — Xq\ > e r . 
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x - e r 


S X "° 


'x +e r x* 



Figure 3. The potential Vq and the intervals I±. 

Here e r > is a small constant we will choose later. We can also compute 

~ 27M 2 w 2 „ 
Aq = , : ,o for a = 0; 



1 



9AM 2 



(4.1) 



y o "(a;o) = -18M 4 (1 - 9AM 2 ) 2 A for a = 0, A = A . 

Letting 

be the principal symbol of P x , we see that po has a nondegenerate hyperbolic critical point 
at (s ,0) and this is the only critical point in the set {po > —ef}. 

4.2. WKB solutions and the outgoing condition. Firstly, we obtain certain approxi- 
mate solutions to the equation P x u = in the region \x — x \ > e r , where V is known to 
be negative. (Compare with |Ra| Sections 2 and 3].) Define the intervals 

1+ = (x + e r , +oo), I = (— oo, x - e r ), I Q = (x - 2e r , x + 2e r ). (4.2) 

Let ipo(x) be a smooth function on I + U I solving the eikonal equation 



sgn(x — Xq 



)v / -^( 



x 



(We will specify a normalization condition for ipo later.) Then we can construct approximate 
WKB solutions 



v£(x; h) = e iMx)/h a^(x; h), u±(x; h) 



-iipa{x)/h 



a ± (x; h), x E I±, 



(4.3) 



{+, — }. Indeed, if 



such that P x u{ = 0{h°°) in C°°(I ± f and a are smooth classical symbols in h, for 7, 5 G 



a 5 7 (x;h) ~J2 hj( 



x 



j>0 



8 Hcnceforth we say that u — 0(h°°) in C°°(I) for some open set /, if for every compact K C I and 
every N, \\u\\c N (i<) — 0(h N ). In particular, this does not provide any information on the growth of u at 
the ends of /. Similarly, we say that u is polynomially bounded in C°°(I) if for every K and TV, there 
exists M such that H^Hc"^) = 0(h~ M ). 
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then the functions have to solve the transport equations 

(2^ (x)d x + <(;r) ± iV x {x))a±W = 



(4.4) 

W Q (x)d x + if%(x) ± ^(x))a^ +1 > = ±i(d 2 x - V 2 (x))a^\ j > 0- 



the latter can be solved inductively in j. We will fix the normalization of later; right 
now, we only require that for x in a compact set, ay ~ 1 in the sense that C~ x < \a^\ < 
C for some /i-independent constant C . Put 

rj = {{x, ±i// (x)) | x G J 7 } C T%, 7 G {+, -}; (4.5) 



then by Proposition 2.1 (with m = 0), 

WF h («?)crJ, 7,5 6 {+,-}. (4.6) 

Now, we show that the Cauchy problem for the equation P x u = is well-posed semiclas- 
sically in I±. For two smooth functions vi,v 2 on some interval, define their semiclassical 
Wronskian by 

W(v 1 ,v 2 ) = v x ■ hd x v 2 - v 2 ■ hd x vi, 

then 

kd x W(v 1 ,v 2 ) = v 2 ■ P x v x - v x ■ P x v 2 . (4.7) 
Also, if W(vi,v 2 ) 7^ and u is some smooth function, then 

W(u, vi)v 2 - W(u, v 2 )vi 
u= ^7? \ • ( 4 - 8 ) 

W(v2,Vi) 

We have W(u±,u±) ~ 1; therefore, the following fact applies: 

Proposition 4.1. Assume that I C M is an interval and U G I is a nonempty open set. 
Let Vi(x; h),v 2 (x; h) G C°°(I) be two polynomially bounded functions such that P x Vj(x; h) = 
0(h°°) inC°°(I) andWivi.v^' 1 is polynomially bounded. (Note that by (4.7), d x W(vi,v 2 ) - 
0(h°°).) Let u(x;h) G C°°(I) be polynomially bounded in C°°(U) and P x u = 0(h°°) in 
C°°(I). Then u = C\V\ + c 2 v 2 + 0(h°°) in C°°(I), where the constants C\, c 2 are polynomi- 
ally bounded. Moreover, Cj = W(u,v 3 ^j)/W(vj,v 3 ^j) +0(h°°). 



Proof. Let Wj = W(u,Vj). Combining (47) and g§, we get \d x Wj\ = 0(/i°°)(|W / i| + |W / 2 |). 
Also, Wj are polynomially bounded on U. By Gronwall's inequality, we see that Wj are 
polynomially bounded on I and constant modulo 0(h°°); it remains to use (4.8). □ 

Now, recall |Dyl[ Section 4] that for X large enough, we have V(x) = V±(e T ±x ) for 
±x > Xq, where A± > are some constants and V±(w) are holomorphic functions in the 
discs {\w\ < e~ A±x °}, and V±(0) = —u>±, where 

u± = (1 + a)((r± + a 2 )(w + ihv) - ah). (4.9) 
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For a and h small enough, we have Reu;± > 0. In |Dyl[ Section 4], we constructed exact 
solutions u±(x) to the equation P x u± = such that 

u± ( x ) = e ±iu±x l h v ± {e^ A±x ) for ± x > X , 

with v±(w) holomorphic in the discs {\w\ < e~ ±x °} and v±(0) = 1. Note that we can use 
a different normalization condition than |Dyl Proposition 4.2], as lmu± = 0(h) under the 
assumptions (1.12). 

Proposition 4.2. For a certain normalization of the functions ipo and 

u± (x) = u+{x) + 0{h°°) in C°°(I±). (4.10) 



In particular, by (4.6) 

WF fc («±| /± ) C T+ (4.11) 



Proof. We will consider the case of u+. By Proposition 4.1, it is enough to show (4.10) for 
±x > Xq, where X Q is large, but fixed. We choose X large enough so that ReV±(w) < 
for | it; | < e~ A±x °. Then there exists a function if)(x) such that 

{d x i>{x)) 2 + V{x) = 0, x > X ; 

ip(x) = oo+x + ip(e~ A+x ), 

with ip holomorphic in {\w\ < e~ A+x °}. We can fix ip by requiring that ^(0) = 0. Take 

u+(x) = e i ^ x)/h a(e- A+x ;h); 

then P x u + = if and only if 

([hA + wD w + A + wi)\w) - uj+} 2 + V)a = 0. 

This can be rewritten as 

-A + {wi)'{w))'a + (2co + + ihA + - 2A + w4>'(w))d w a + ihA + wd 2 w a = 0. 

We will solve this equation by a power series in w and estimate the terms of this series 
uniformly in h. Let us write 

^'(w;h) = ^^ipi(h)w l , a(w;h) = Tj aj{h)w^ 

l>0 j>0 

and solve for a, with the initial condition qq = 1, obtaining 

a ^ h) = U + l)(2u + i\hA + (j + l)) E Mh)(l + 2j-l)a^(h). 

We claim that for some R, all j, and small h, \aj(h)\ < Pj . Indeed, we have \2u + +ihA + (j + 
1)| > e > 0; combining this with an estimate on ipi, we get 

\a j+1 \ < S\l + 2j - l)\a^\ <2CJ2 ^K'-d 

o</<i o<i<j 
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for some constants C and S. We can then conclude by induction if R > 2C + S. In a similar 
way, we can estimate the derivatives of aj in h; therefore, a(w; h) is a classical symbol for 
\w\ < Br 1 . 

Now, we take Xq large enough so that e~ A+x ° < R^ 1 and restrict ourselves to real 
x > Xq. We can normalize ?/>o so that ip(x) = ipo(x) + hipi(x; h) for some classical symbol 
ipi ; then 

u+(x) = jM^/h^M^a^-A+x. 
The expression in square brackets is a classical symbol; therefore, this expression solves the 



transport equations (4.4); it is then equal to a constant times a+, modulo 0(h°°) errors. □ 



4.3. Transmission through the barrier. First of all, we establish a microlocal normal 
form for P x near the potential maximum. Let Eq > be small; define 

K = {\x-x \ <e , |f | <e }cT*R. 
We pick e r small enough, depending on e , such that e r < e /2 and 

{ Po = o}cK u{Jr s v 

7,5 



4.1 



that e r controls how close we are to the 



with defined in (4.5). (Recall from Section 
trapping region.) We also assume that e is small enough so that (x , 0) is the only critical 
point of p in K . 

Proposition 4.3. For Eq small enough ande r small enough depending on Eq, there exists a 
symplectomorphism $ from a neighborhood of Kq onto a neighborhood of the origin in T*IR 



and operators B±, B 2 quantizing $ near K x $(i£o) in the sense of Proposition 2.5, such 
that P x is intertwined via (Bi, B 2 ) with the operator SQ(f3) microlocally near Kq x <&(Kq), 
with S E elliptic in the class *°(R) ; 



Q(P) = hxD x - 0, 
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and (3 = (3(u, u, A, jl, k; h) a classical symbol. Moreover, the principal part (3q of (3 is real- 
valued, independent ofu,jl, and vanishes if and only A = \ (u,k). Also, 

V (x ) 



Pa 



Po = - , 

y/-2Vf(x ) 

27Mfa 2 - A(l - 9AM 2 ) 



+ O(V (x ) 2 ), 
+ O(V (x ) 2 ) fora = 0. 



(4.12) 
(4.13) 



2V A(l - 9AM 2 ) 

Finally, $(K ) D K Q = {\x\ < i , |£| < e } for some e depending on e , and, with I 
defined in (4.2), 

$(r+) c f + = {±x > 0, £ = Po/x, |f | < eo/2}, (4.14) 

$(r ± ) cfi = {T(>0,x = /V£, |x| < £ /2}; (4.15) 

{p = 0} fl {s G Jo} C ^-^^o) C tf . (4.16) 



Proof. First of all, we use |CdVPat Theorem 12] to construct &,Bi,B 2 conjugating P x 
microlocally near the critical point (xq, 0) to an operator of the form f[hxD x ], for some 



symbol f(s;h), where the latter employs the formal functional calculus of Section 2.2 



The techniques in the proof are similar to those of Section 3^ of the present paper, with 
appropriate replacements for Propositions |2. 10] and 3.7 therefore, the proof goes through 
for complex valued symbols with real principal part. 

Let f be the principal part of /; then p o = fo(x£). Note that $(x ,0) = (0,0). 
The level set {p = V (x )} at the trapped energy contains in particular the outgoing 
trajectory {x > x , £ = a/Vo(x ) — V (x)}; we can choose $ mapping this trajectory into 
{x > 0, £ = 0}. Since the latter is also outgoing for the Hamiltonian flow of x£, we have 
d s fo(0) > 0; it follows that d s fo(s) > for all s (if d s fo vanishes, then p has a critical 
point other than (xq,0)). The function f(s;h) not uniquely defined; however, its Taylor 
decomposition at s = 0, h = is and we can compute in particular 

fo(s) = V (x ) + s^/-2V£(x ) + 0(s 2 ). (4.17) 

Therefore, for e r small enough, we can solve the equation f(s; h) = for s; let (3 be 
the solution. We now write f(s;h) = fi{s;h)(s — (3) for some nonvanishing fi and get 
f[hxD x ] = SQ((3) microlocally in <&(K ), with S = fi[hxD x ] in <&(K ) and extended to be 
globally elliptic outside of this set. The equation (4.12) follows from (4.17), while (4.13) 



follows from ( |4.12[ ) and ( |4.1| ). 

Finally, po = on each and thus x£ = (3 on $(T*). By analysing the properties of $ 
near (xq, 0), we can deduce which part of the sets {x£ 
follow for e r small enough. 



/3}eachr* maps into; (4.14)-(4.16) 



□ 



We now describe the radial quantization condition and provide a non-rigorous explana- 
tion for it. Recall from |Dyl[ Section 4] that (u, A, k) is a pole of R r if and only if the 
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functions u±, studied in the previous subsection, are multiples of each other. Assume that 
this is true and u = u + ~ it_. However, by (4.11) the function B\u is microlocalized on the 
union of Tjjj, but away from r±; it also solves Q(/3)Biu = microlocally. By propagation 
of singularities, this can happen only if the characteristic set of Q(f3) is {x£ = 0}, and in 
this case, B\u is smooth near x = 0. Then B\u must be given by x J/3///l , with (3 G — ihN and 
N denoting the set of nonnegative integers. Therefore, we define the radial quantization 
symbol J-" r (m, u, u, k\ h) as the solution A + ihjl to the equation 

/3(d), z>, A, fl, k] h) = —ihm, m G Z, < m < C m . 



The expansions for T r near a = described in Proposition 1.5 follow from (4.13) and (B.12). 

We now prove the rest of Proposition 1.5 We start with quantifying the statement that 
in order for the equation Q(/3)u = to have a nontrivial solution smooth near x = 0, the 
quantization condition must be satisfied: 

Proposition 4.4. Assume that (5 e C satisfies 

|/?| < Cp, I Im/3| < Ceh, min |/3 + ihm\ > Cfh, (4.18) 

for some constant . Let U C R 6e a bounded open interval, I C U be a compact interval 
centered at zero, and X e \l/^° mp (]R). T/ien i/iere exzst constants C and N such that for 
each u G L 2 (R), 

\\Xu\\ L 2 (I) < Ch- N \\Q((3)Xu\\ LHu) + 0(/i°°)||«|| i2(R) . (4.19) 

Proof. First, assume that Im/3 > h. Let /' be an interval compactly contained in / and 
centered at zero. We will use the fact that every C°°(I') seminorm of Xu is bounded by 
Ch~ N \\Xu\\i,2(i) + 0{h co )\\u\\i J 2^ for some constants C and N, depending on the seminorm 
chosen. (Henceforth C and iV will be constants whose actual values may depend on the 
context.) Since Xu G C°°, we can write 

hD x {x~ il3/h Xu) = x~ x - ifilh Q($)Xu. 

However, Re(— i/3/h) > 1; therefore, x~ % $l h Xu vanishes at x = and we can integrate to 
get 

\\x-W h Xu\\ L ~w < C\\Q{/3)Xu\\ L 2 {I) . (4.20) 

On the other hand, 

\\Xu\\ Lac{n < Ch- No \\Xu\\ L 2 {I) + 0(h°°)\\u\\ L 2 m . (4.21) 



for some constants C and A^ . Now, take a large constant x; using (4.21) in {\x\ < h*} 
and (4.20) elsewhere, we get 

\\xu\\ LHI) < + ^-^HQ^xwii^^) + o^ 00 )!^!!^^); 



taking x large enough, we get (4.19). 
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For the general case, we choose an integer M large enough so that Im(/3 + ihM) > h. 
Let v be M-th Taylor polynomial of Xu at zero. Since 

\\Q{P + ihM)D A JXu\\ LHl) = \\D^Q{P)Xu\\ L 2 {I) < Ch- N \\Q(J3)Xu\\v {u) + 0(h°°)\\u\\v m , 

we apply the current proposition for the case Im /3 > h considered above to get 

\\Xu-v\\ LHl) < C\\D^Xu\\ LHl) < Ch- N \\Q(f3)Xu\\ L2(u) + 0(h^)\\u\\ L , m] 

therefore, ||Q(/3)u is bounded by the same expression. However, one can verify directly 
that if (3 is C^ l h away from —ihN, then 

IMb(i) < Ch- N \\Q((3)v\\ L 2 {I) - 
this completes the proof. □ 

Now, we show that each radial pole lies within o(h) of a pseudopole: 



Proposition 4.5. Assume that /3(h) satisfies (4.18). Then for h small enough, (uj,X } k) is 
not a radial pole, and for each compact interval /cK, there exist constants C and N such 
that 

|| ljR r (co, A, k)lf\\i'2^ L 2 < Ch~ N . 

Proof. Let u E H 2 oc (M.) be an outgoing solution to the equation P x u = f, with / G L 2 
supported in a fixed compact subset inside the open interval I. Then 

u(x) = c±u±(x), ±x ^> 0, 

for some constants c±. Clearly, \c±\ < C\\u\\i?m- Using the method of proof of Proposi- 
tion 44, we get 

\\u - c ± u ± \\ L 2 (InI±) < ChT l \\f\\ L ,. (4.22) 

Next, let x £ Co°(^o) be equal to 1 near the complement of I + U I and Bi be the operator 
introduced in Proposition 4.3 consider the compactly microlocalized operator 

T = SQ{fi)B x - B X P X . 



Then by Proposition 2.5, we can write T = TX + 0(h°°), where X 6 v]/™ mp is a certain 
operator vanishing microlocally on Kq. By (4.16), we can further write X = X1+X2, where 
Xj G ^- m P, WF h {Xj) n K = 0, WFhiXi) n {po = 0} = and WF h (X 2 ) C {x e 1+ U !_}. 
By ellipticity, 

ll^i«IU 2 <c||/|| L2 + o(OIHI^(/)- 

Take X ± € Cg°(I±) such that x± = 1 near J ± n K(WF h (X 2 )) (here vr : T*IR -> R is the 
projection map onto the base variable and 7r(WF/ l (X 2 )) is a compact subset of 1+ U 
then by fl4~22j ), 

||X 2 (u - Ml )|| L2 < CTT^/H^ + 0(^)|| M || L2(/) , 
ui = c + x+(x)u + + c_x-(a:)u_. 
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It follows that 

\\SQ(^B 1 u-TX 2 u 1 \\ L ,<Ch- 1 \\f\\ L2 + 0(h^)\\u\\ L 2 {r) . 
Combining Q4.11D with Q4.14D and the fact that WF h (X 2 ) nK = (ft, we get 

WF h (TX 2X ±(x)u ± ) GT+\K . 

The projections of the latter sets onto the x variable do not intersect I Q = {\x\ < e }; 
therefore, for some open Uq containing Iq, 

\\TX 2Ul \\ L2{Do) = 0^)11^1^. 

Using ellipticity of S, we then get 

\\Q(P)B lU \\ L2{Uo) < Ch^WfWv +0(h°°)\\u\\ LHl) . 



Applying Proposition |4.4| to B\u on J and using that B\ is compactly microlocalized, we 

get 

\\XB 1 u\\ L 2<Ch- N \\f\\ I? +0{h^)\\u\\ L 2 {I) , 



for any X G \l/™ mp microlocalized in Kq. Using the elliptic estimate and (4.16), we get 

\\u\\ L 2 {Io) <Ch- N \\f\\ L i + 0(h™)\\u\\ L i {I) . 



From here by (4.22), 



\c ± \ < C||c ±U± || L2(J±n/o) < C7i- JV ||/|| L2 + 0(h°°)\\u\\ L 2 (I y, 



combining the last two estimates with (4.22), we get the required estimate: 

IMU/)< ^-"11/11^ + o(h°°)\\u\\ L 2 {I) . □ 



To finish the proof of Proposition 1.5, it remains to show 



Proposition 4.6. Fix uj,v,k satisfying (1.12), m G N bounded by a large constant C m , 
and let V be the set of all X such that 

\/3(uj, v, k, A, /i; h) + ihm\ < h/3. 

Then for h small enough, R r (u,X,k) has a unique pole Ao in V, and Xq is within 0(h°°) 
of T r {rn,uj, u,k;h). Moreover, we can write 



R r (u, A, k) 



where the family of operators S(X) : L 
S(X ) is a rank one operator. 



comp v 



S(X) 

A-A ' ' l - " 
-+ Lie 



XeV, 

(M) is bounded polynomially in h and 
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Proof. We will use Proposition 4.3 and the fact that = 0(h) to extend the WKB solutions 



from I± to the whole 



Consider the locally integrable functions 

(x ± iO) ifS/h . 



u 



± 



solving the equation Q(j3)u± 
basic properties of (x ± i0) b .) 



= 0. (See for example |Ho[ Section 3.2] for the definition and 
We have 



WF h (u+) C {£ = 0} U {x 



U ± {X) 



X 



iP/h 



microlocally near {x > 0, £ = 0}, 



= 0, ±£>0}; 

( x ) = e T*P/h(_ 



X 



AP/h 



microlocally near 



{x < 0, £ = 0}. Using the formulas for the Fourier transform of u± |Hol Example 7.1.17], 

we get 



j l iP/h e TP7r/(2h) 



U±{X) 



-1-iP/h ±ix£/h 



±£ G K^}, for every x C^°(0, oojsuch that x 



microlocally near {x = 

.fQ C (0, oo). Let B 2 be the operator constructed in Proposition 
0(fr°°) in C°°(Iq), and B 2 u\ = c ± u^ + 0(h°°) in l7°°(7 ± n J ) forsbme constants c± 



4.3 



1 near 

By (|4T6l) P XJ B 2 ui = 
1. 



To prove the latter, we can use the theory of Fourier integral operators and Lagrangian 
distributions to represent B 2 u± in the form ( 4.3[ ) microlocally near rjjj; the symbols in these 
WKB expressions will have to solve the transport equations. Then we can use B 2 u± to 
extend u± to I± U Jq so that P x u± = 0(h°°) there. We claim that 



u 



± 



c ±l u% + c ±2 T(-if3/hy l u~ + 0(h c 



(4.23) 



in C°°(I T PI Jo); with c±j constants such that cj-j and c ± j are polynomially bounded in h. 



To show (4.23), we can apply the theory of Lagrangian distributions to B 2 u± one more 
time; alternatively, we know that this function is an 0(h°°) approximate solution to the 
equation P x u 

and T^. 
with P x u% = 0(h°°) in C 
using either of the families 
get u± = ut + 0(h°°) in C° 



on Iq fl J=p, and we have control on its L 2 norm when microlocalized to 
Thus, we can extend u± to the whole R as a polynomially bounded family 



and ( |4.23 ) holding on I T . Similarly we can extend u 



u±,u±) 



It now follows from (4.23) that 



in Proposition 4.1 together with Proposition 4.2 



we 



W(u+,u-) = c(T(- i f3/h)- 1 + 0(h°°)), 
with c and c -1 bounded polynomially in h. By |Dyl[ (4.8)], we get 



R r (u, A, k) 



S(cj, A, k) 
W{u 



u. 



with the family S holomorphic and bounded polynomially in h. Moreover, for W(u + , uJ) = 
0, S is proportional to u + ® u + and thus has rank one. We are now done if we let A be 
the unique solution to the equation W(u + ,u~) = in V. □ 
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Appendix A. Grushin problems for several commuting operators 

A.l. Global Grushin problem. Assume that P\, . . . ,P n are pseudodifferential operators 
on a compact manifold M, with Pj G \I/ fcj (M) and kj > 0. 



Definition A.l. We say that X = (Ai,...,A n ) G C n belongs to the joint spectrurr^j of 
Pi, ... , P n , if the joint eigenspace 

{u G C°°(M) | P jU = XjU, j = 1, . . . , n} 

is nontrivial. (In our situation, one of the operators Pj will be elliptic outside of a compact 
set, so all joint eigenf unctions will be smooth.) 



The goal of this appendix is to extract information about the joint spectrum of Pi, . . . , P n 
from certain microlocal information. Essentially, we will construct exact joint eigenfunc- 
tions based on approximate eigenfunctions and certain invertibility conditions. The latter 
will be given in the form of operators A±, . . . , A n , with the following properties: 

(Gl) Each Aj can be represented as A 1 - + A", where A'- is compactly microlocalized and 

has operator norm 0(h~ r ); A" G h~ r ^f~ k3 (M). Here r > is a constant. 
(G2) The commutator of any two of the operators Pi, ... , P n , Ai, . . . , A n lies in h°°^f~°°(M). 

We would like to describe the joint spectrum of Pi, . . . , P n in a ball of radius o{h r ) centered 
at zero. First, we consider a situation when there is no joint spectrum: 

Proposition A.l. Assume that conditions (Gl) and (G2) hold and additionally, 

n 

A i p j = 1 mod h°°^f~ 00 (M). 

Then there exists 5 > such that for h small enough, the ball of radius 5h r centered at zero 
contains no joint eigenvalues of Pi, . . . , P n . 

Proof. Assume that u G L 2 (M) and PjU = XjU, where |Aj| < 5h r . Then 

n n 

= J2 M P 3 - X i) U = ( J + h 00 ^- 00 ^ - X 3 A 3 U - 
3=1 3=1 

It follows from condition (Gl) that ||y4j||^2_ s>L 2 = 0(h~ r ); therefore, 

\\u\\ L 2 = 0(5 + h°°)\\u\\ L 2 
and we must have u = for 5 and h small enough. □ 



9 Strictly speaking, this is the definition of the joint point spectrum. However, the operators we study 
in Section [3] are joint elliptic near the fiber infinity, as in Proposition A. 4 thus all joint spectrum is given 
by eigenvalues. 
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Now, we study the case when the joint spectrum is nonempty. Assume that S x : C — > 
C°°(M) and S 2 '■ T>'(M) — > C are operators with the following properties: 

(G3) Each Sj is compactly microlocalized with operator norm 0(1). 
(G4) S 2 S 1 = l + 0{h°°). 

(G5) If Q is any of the operators P x , . . . , P n , A x , . . . , A n , then QSi G /i 00 ^- 00 and S 2 Q G 
(G6) We have 

n 

A j p j = 1 - SiS 2 mod h^-^iM). 

i=i 

Note that (G5) implies that the image of Si consists of 0(/i oc, )-approximate joint eigen- 
functions. For n — 1, one recovers existence of exact eigenfunctions from approximate 
ones using Grushin problems, based on Schur complement formula; see for example |HiSj 



Section 6]. The proposition below constructs an analogue of these Grushin problems for 
the case of several operators. This construction is more involved, since we need to combine 
the fact that Pj commute exactly, needed for the existence of joint spectrum, with microlo- 
cal assumptions (G1)-(G6) having 0(h°°) error. Note also that condition (G2) does not 
appear in the case n = 1. 

Proposition A. 2. Assume that the conditions ( Gl )-( G6) hold and the operators Pi, . . . , P n 
commute exactly; that is, [Pj,Pk] = for all j,k. Then there exists 5 > such that for 
h small enough, the ball of radius 5h r contains exactly one joint eigenvalue of Pi, ... ,P n . 
Moreover, this eigenvalue is 0(h°°) and the corresponding eigenspace is one dimensional. 

Proof. We prove the proposition in the case n = 2 (which is the case we will need in the 



For A 



T(A) 



V 



Si 






: Hi — ^ H2 



present paper); the proof in the general case can be found in Appendix A. 3 
(Ai, A 2 ) G C 2 , consider the operator 

fPi-Xi -A 2 
P 2 - A 2 Ai 
S 2 

s 2 

Ui = L 2 {M) © H h kl - k2 {M) © C 2 , H 2 

The conditions (G1)-(G6) imply that for 

/ Ai A 2 
-P2 Pi 
S 2 

s 2 



Q 



\ 



Si 






Si 




o\ 

Si 


°/ 

H h kl (M) @H~ k2 (M) ©C 2 . 



: 1~L 2 — y Hi, 



we have T(0)Q = I + O n ^ m {h°°) , QT(0) = J + O^^/i 00 ). By (Gl) and (G3), we have 
||Q||« 2 ->'Hi = 0(h~ r ); therefore, if 5 > and h are small enough and |A| < 5h r , then T(A) 
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is invertible and 



Now, let |A| < 5h r and put 



HT(A)- 1 ! 



0{h~ r ). 



(u(X),u 2 (X),f(X)) = T(A) (0, 0, 1, 0), /(A) = (/i(A) J 2 (A)) e C 2 . 

This is the only solution to the following system of equations, which we call global Grushin 
problem: 



(Pi - Ai)u(A) - A 2 u 2 {X) + 5i/i(A) = 0, 
(P 2 - A 2 )u(A) + A lU2 {X) + SJ 2 (X) = 0, 

S 2 u(X) = 1, 
S 2 u 2 (X) = 0. 



(A.l) 
(A.2) 
(A.3) 
(A.4) 



We claim that A is an element of the joint spectrum if and only if /(A) = 0, and in that 
case, the joint eigenspace is one dimensional and spanned by u(X). First, assume that u is 
a joint eigenfunction with the eigenvalue A. Then T(X)(u, 0, 0, 0) = (0,0, s, 0), where s is 
some nonzero number; it follows immediately that /(A) = and u is a multiple of u(X). 

Now, assume that /(A) = 0; we need to prove that u(X) is a joint eigenfunction for 
the eigenvalue A. By ( |A.l ) and (A.2), it suffices to show that u 2 (X) = 0. For that, we 
multiply (A.2) by Pi — Ai and subtract (A.l) multiplied by P 2 — A 2 ; since /(A) = and 
[Pi,P 2 ] =0, we get 

((Pi - Xi)A x + (P 2 - X 2 )A 2 )u 2 (X) = 0. 

Recalling (G6), we get 

(/ - SiS 2 + O n - kl - k , H - kl - k2 (5 + h°°))u 2 (X) = 0. 

h h 

By flAl| , (/ + 0(5 + h°°))u 2 (X) = and thus u 2 (X) = 0. The claim is proven. 

It remains to show that the equation /(A) = has exactly one root in the disc of radius 
Sh r centered at zero, and this root is 0(h°°). For that, let QT(X) = I — R{X); we have 

/ Ai^i + X 2 A 2 0\ 

-AiP 2 + A 2 Pi 

AiS 2 

A 2 S 2 0/ 

T(A)' 1 = (J + P(A) + (J - P(A))- 1 P(A) 2 )g. 
One can verify that P(A) 2 Q(0, 0, 1, 0) = Hl (h°°) and then 

f(X) = X-g(X;h), 



R(X) 



\ 



+ 



where g(X; h) = 0(h c 
principle. 



uniformly in A. It remains to apply the contraction mapping 

□ 
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Finally, we establish a connection between global Grushin problem and meromorphic 
resolvent expansions, using some more information about our particular application: 

Proposition A. 3. Assume that n = 2, Pi, P 2 satisfy the properties stated in the beginning 
of this subsection, [Pi,P 2 ] — 0, ki > 0, and P\ — X is elliptic in the class for some 
A G C. If V is the kernel of Pi, then by analytic Fredholm theory (see for example |EvZw[ 
Theorem D.4]j, the resolvent 

R(\) = (P 1 - A) _1 | y : H~ kl (M) n V -> L 2 (M) n V 

is a meromorphic family of operators in A G C with poles of finite rank. Then: 



(1) Assume that the conditions of Proposition A.l hold and let 5 > be given by this 
proposition. Then for h small enough, R(X) is holomorphic in {|A| < 5h r } and 
\\R(\)\\ L 2 nV ^ L 2 = 0{h~ r ) in this region. 



(2) Assume that the conditions of Proposition A. 2 hold and let (Ao, A 2 ) be the joint 
eigenvalue and 5 > the constant given by this proposition. Suppose that A° = 0. 
Then for h small enough, 



R(X) = S(X) + 



n 



A-A ' 



IAI < 6h r 



where S(X) is holomorphic, U is a rank one operator, and the L 2 fl V — > L 2 
of 5(A) and U are 0(h~ N ) for some constant N . 



norms 



Proof. 1. We have 

A 1 (P 1 -X)=I-XA 1 + h 00 ^- 00 (M)onV; 
the right-hand side is invertible for 8 small enough. Therefore, R has norm 0(h~ r ). 

2. We know that R(X) has a pole at A if and only if there exists nonzero u G L 2 (M) fl V 
such that (Pi — X)u = 0; that is, a joint eigenfunction of (Pi,P 2 ) with joint eigenvalue 
(A,0). Therefore, Ao is the only pole of R(X) in {|A| < 5h r }. 

Now, take A 7^ A , |A| < 5h r , and assume that v G H h hl (M) fl V and u = R(X )v G 
L 2 (M) D V. Let T(A) be the family of operators introduced in the proof of Proposition A. 2 
with Ai = A and A 2 = 0; we know that T(A) is invertible. We represent T(A) -1 as a 4 x 4 
operator- valued matrix; let (A) be its entries. We have T(X)(u, 0, 0, 0) = (v, 0, c, 0) for 
some number c. However, then (w, 0, 0, 0) = T(A) -1 (t> , 0, c, 0); taking the third entry of this 
equality, we get T 3 "[ 1 (A)u-|-r3^ 1 (A)c = 0. Now, T 3 J 5 1 (A) = /i(A), with the latter introduced in 



the proof of Proposition |A.2| Therefore, we can compute c in terms of u; substituting this 
into the expression for u, we get the following version of the Schur complement formula: 



P(A)= Tf/(A) 



/i(A) 



(A.5) 



Next, by the proof of Proposition A. 2, fi(Xo) = and /i(A) = A + 0(h°°). Therefore, we 
may write /i(A) = (A — Xo)/g(X), with g holomorphic and bounded by 0(1). Let uq be 
the joint eigenfunction of (Pi,P 2 ) with eigenvalue (Aq,0); then LT = — 5f(Ao)T 1 3 1 (Ao)T 3 ^ 1 (Ao) 
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is a rank one operator, as T 13 1 (A ) acts C — > V and Hu ■ 
operators T^ 1 are polynomially bounded in h, we are done. 



(1 + O(h°°))u . Since the 

□ 



A. 2. Local Grushin problem. In this subsection, we show how to obtain information 
about the joint spectrum of two operators Pi, P 2 based only on their behavior microlocally 
near the set where neither of them is elliptic. For that, we use global Grushin problems 
discussed in the previous subsection. Assume that Pi G ^^(M), P 2 G ^f^(M) satisfy 

(El) The principal symbol pj of Pj is real-valued. 

(E2) The symbol p w is elliptic in the class S kl (M) outside of some compact set. As a 
corollary, the set 



is compact. 

Next, assume that A±,A 2 are compactly microlocalized operators on M such that: 

(LI) For each j and every bounded neighborhood U of K, Aj can be represented as 
A'j + Aj, where both A'j and A" are compactly microlocalized, ||^|| = 0(h~ r ), 
WF h (A'j) C U x U, and A'j G h~ r ^l° mp (M). Here r > is some constant. 

(L2) The commutator of any two of the operators Pi, P 2 , A 1: A 2 lies in h co ^~ co (M). 

Finally, let Si : C — > C°°(M), S 2 : V'(M) — > C be compactly microlocalized operators such 



(L3) || Sj\\ = 0(1) and WF h (Sj) C K. 
(L4) S 2 Si = l + 0(h°°). 

(L5) If Q is any of the operators P u P 2 , A u A 2 , then QS 1 G h 00 ^" 00 and S 2 Q G h 00 ^- 00 . 



Proposition A. 4. 1. If the conditions (E1)-(E2) and (L1)-(L2) hold, and 

A X P X + A 2 P 2 = I 

microlocally near K x K , then there exists 5 > such that for h small enough, there are 
no joint eigenvalues of P\,P 2 in the ball of radius 5h r centered at zero. 

2. If the conditions (E1)-(E2) and (L1)-(L5) hold, [Pi,P 2 ] = 0, and 



microlocally near K x K , then there exists 5 > such that for h small enough, the ball 
of radius 5h r centered at zero contains exactly one joint eigenvalue A of P\,P 2 . Moreover, 
A = 0(h°°) and the corresponding joint eigenspace is one dimensional. 

Proof. We will prove part 2; part 1 is handled similarly. Take small e > and let Xe £ 
Cq > (M) be supported in (—e,e) and equal to 1 on [— e/2, e/2]. Also, let t/j £ G C°°(R) satisfy 
tip £ (t) = 1 — Xe(t) f° r & U then ip e {t) = for \t\ < e/2. The function ip e is a symbol of 
order —1, as it is equal to t^ 1 for \t\ > e. 



K = {(x, G T*M I p w (x, = p 20 (x, = 0} 



that: 



A 1 P 1 + A 2 P 2 = / - 



(A.6) 
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By (El), we can define the operat ors x JP,-], ipe\Pj] G W™(T*M) usin g the formal func- 



tional calculus introduced in Section 



2.2 



By (E2) and Proposition |2.4|^ e [Pi] G ^ d kl (M), 
and Xe[-Pi] G ^° mp . Therefore, we can define uniquely up to /i 00 * -00 the operators 

Xe = Xe[Pl]Xe[P2] G ^(M), *I>e[ P l] G ^ W. X-frhM^] G ^(M). (A.7) 

By Proposition |2.3[ these operators commute with each other and with P\,P 2 modulo 



h°°ty °°. Let y be any of the operators in (A.7); we will show that it commutes with each 



Aj modulo h°°ty 00 . Take a neighborhood U of K so small that \pw\ + \p2o\ < e/4 on 
£/; then y is either zero or the identity operator microlocally on U. By (LI), decompose 
A 3 - = A'j + A], where WF h (AJ) C E7 x £/ and ^ G *™ mp . We have Aj = A'j microlocally 
away from U x U; therefore, [A'', Pfc] = microlocally near T*M \ U . By Proposition 



2.3 



[Aj,y] = microlocally near T*M \ U; therefore, the commutator [A,-,y] is compactly 
micro localized and WF^QA,-, Y]) C U X U. However, since Y = or Y = I microlocally in 
U, we have [Aj,Y] G h 00 ^' 00 , as needed. 

Since X £ = I microlocally near K and WFh(Sj) C K, we get (/ — X s )Sx, S 2 (I — X e ) G 
/i°°\I/ _00 . Multiplying ( |A.6 ) by X £ , we get for e small enough, 



(X £ A 1 )P 1 + (X £ A 2 )P 2 + S 1 S 2 = X £ mod/1 00 *- 00 . (A.8) 



Next, by Proposition 2.3 



A[Pi}Pi + Xe[Pi}A[P2}P2 = I~X £ mod/i 00 *" 00 . (A.9) 
Adding these up, we get 

{X e A X + MPl\)Pl + (XeA 2 + Xe[Pl]lpe[P2])P2 + S X S 2 =1 + h™^ 00 . (A.10) 

The operators P x , P 2 , A x = X E A X + ^ e {P x ], A 2 = X £ A 2 + Xe{Pi]^e[Pi], S u S 2 satisfy the 



assumptions of Proposition A. 2 Applying it, we get the desired spectral result. □ 



A. 3. Proof of Proposition A. 2 in the general case. In this subsection, we prove 



Proposition A. 2 for the general case of n > 2 operators. For simplicity, we assume that 



k\ — • ■ • — k n — 0; that is, each Pj lies in * (M). (If this is not the case, one needs to 
replace L 2 (M) below with certain semiclassical Sobolev spaces.) 

Let V be the space of all exterior forms on C n ; we can represent it as Vsven © Vodd, where 

VWen = A 2J C\ Vodd = A 2i+1 C n 

j>0 j>0 

are the vector spaces of the even and odd degree forms, respectively. Note that Vsven and 
Vodd have the same dimension. Define the spaces 

L| ven = L 2 (M) ® V Even , L 2 Qdd = L 2 (M) <g> Vbdd, L 2 V = L 2 (M) <g> V. 

We call elements of L v forms. They posess properties similar to those of differential forms; 
beware though that they are not differential forms in our case. We will use the families of 
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operators (Aj) and (Pj) to define the operators 
given by the formulas 



d P ,d* A :L v ^L 2 



Vi 



d P (u ®v) = ^(Pju) <g> (ej A v), 



n 



d* A (u <g) v) = 2^(A,-w) g) (i ei u); 
w G L 2 (M), D6K 

Here ei, . . . , e n is the canonical basis of C n . The notation i e . is used for the interior product 
by ej] this is the adjoint of the operator v h-> ej A v with respect to the inner product on 
V induced by the canonical bilinear inner product on C n . Note that dp and d* A map even 
forms to odd and vice versa. 

A direct calculation shows that under the assumptions (G1)-(G6), 

(dp + d* A f = 1- SA ® I v + Oy-^h 00 ). (A.ll) 

Here Jy is the identity operator on V, while / is the identity operator on L v . Moreover, 
since the operators Pi, . . . , P n commute exactly, we have 

4 = 0. (A. 12) 

For A = (Ai, . . . , A ra ) G C n , define the operator 

{ ' V S * ® I V ) ■ Ul ^ 

"Hi = -^Evon © ^Odd, %2 = -^odd © ^Even- 
Here <ip_A is defined using the operators Pi — Ai, . . . , P n — X n in place of Pi, . . . , P n . It 
follows from (A.ll) that for 

Q= ((d P + d* A )\ Lhdd S X ®I V \ 

we have QT(0) = I + O ni ^ ni (h°°), T(0)Q = 1 + H2 ^ H2 {h°°). Moreover, it follows from 
(Gl) and (G3) that \\Q\\ H ^ ni = 0{h~ r ). Therefore, for |A| < 5h r and h and 5 > small 
enough, the operator T(X) is invertible, with ||T(A) _1 ||^ 2 ^ 1 = 0{h~ r ). 

Assume that |A| < 5h r and let 1 G Vsven be the basic zero-form on C n . Put (ot(A), v (A)) = 
T(A) _1 (0, 1), where a(A) G £|ven5 V W e K)ddj then (a(A), u(A)) is the unique solution to 
the system 

(rf P _ A + <£)a(A) + 5i(l) ® «(A) = 0, 

(5 2 <8> Jy)a(A) = 1. ( ' ' 
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We further write v(X) = f(X) + w(X), where /(A) is a 1-form and w(X) is a sum of forms of 
degree 3 or more. Note that both / and w are holomorphic functions of A, with /(A) G C n . 

We claim that A is a joint eigenvalue of Pi, . . . ,P n if and only if /(A) = 0. First of all, 
if u is a joint eigenfunction, then T(X){u ® 1, 0) = c(0, 1) for some scalar c ^ 0; therefore, 
/(A) = and the joint eigenspace is one dimensional. 

Now, assume that /(A) = 0. We will prove that the solution to (A. 13) satisfies a(X) = 
u (g) 1 for some u G L 2 (M); it follows immediately that (Pi — X\)u = ■ ■ ■ = (P n 
Let a = u £g> 1 + /3, where /3 is a sum of forms of degree 2 or higher. Then by (A. 12), 



(dp_ x + d* A ) 2 {u <g> 1) G L 2 (M) (g) 1. 



A n )w = 0. 
"1, 

(A.14) 



Next, we get from (A. 13) 



(S 2 ® Iv)(dp-A + d* A )a + (1 + 0(/i°°))v = 0. 

The components of this equation corresponding to odd forms of degree 3 or higher depend 
only on (3 and w; therefore, for h small enough, w = W(3 for some operator W of norm 
0{h~ r ). Since / = 0, we get v = W{3; therefore, by (A.14) and (A. 13) multiplied by 
dp-x + d* A , 

(rf P _ A + d\ff3 + (rf P _ A + d A ){S x {l) ® W(3) G L 2 (M) ® 1. 

Taking the components of this equation corresponding to forms of even degree 2 or higher 
and recalling (A.ll), we get 

((/ - S^i) ®I V + 0(5 + h°°))f3 = 0. 



However, (S 2 ® I V )P = by ( |A.13p ; therefore, 

(/ + 0(5 + /i°°))/3 



0. 



It follows that /3 = and the claim is proven. 

It remains to show that the equation /(A) = has exactly one solution in the disk of 
radius 5h r . For that, we write QT(X) = I — R(X), 

T(A)- 1 = (I + R(X) + i?(A))- 1 i?(A) 2 )Q. 

We have Q(A)(0, 1) = (^(1) ® 1, 0) and 



R(X) 



(d P + d* A )d x 
(S 2 ®I v )d x 



0(/i c 



Here d\ is constructed using Ai, . . . , A n in place of Pi, ... , P n . Now, we use that P(A) 2 Q(0, 1) 
Wl (/i°°) to conclude that f(X) — X — g(X;h) with g = 0(h°°); it then remains to use the 
contraction mapping principle. 
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Appendix B. Numerical results 

B.l. Overview. This section describes a procedure for computing the quantization symbol 
J-"(m, /, k) from Theorem [I] to an arbitrarily large order in the case 

I' = I - \k\ = 0(1). (B.l) 

The reason for the restriction I' = 0(1) is because then we can use bottom of the well 
asymptotics for eigenvalues of the angular operator; otherwise, we would have to deal 
with nondegenerate trajectories, quantization conditions for which are harder to compute 
numerically; see for example [CdVj. 



We first use the equation (1.16); once we get rid of the semiclassical parameter h (re- 
membering that the original problem was /i-independent), the number u = J-"(m, l,k) is 
the solution to the equation 

g r (m,u,k) = g e (l',uj J k). (B.2) 

Here Q r ,Q e are the non-semiclassical analogues of J-" r , T e \ namely, (1.13) and ( 1.15[ ) take 
the form 

A = Q r (m, u, k) ~ Gj(m, u, k), 
j>o 

\ = g 9 (i',u,k)~J2G°(i'^ik) 1 

respectively. The functions Qj, are homogeneous of degree 2 — j in the following sense: 
£}(m, M s u, sk) = s^Q^m, to, k), Q]{1\ M s u, sk) = s^Q^l', u,k), s > 0. (B.3) 

Here M s u = sKeu + i Im u; the lack of dilation in the imaginary part of u> reflects the fact 
that it is very close to the real axis. 



We will describe how to compute for an arbitrary value of j in Section B.3 

The method is based on a quantization condition for barrier-top resonances, studied in 



Section 4.3 their computation is explained in Section B.2 and a MATLAB implementation 



and data files for several first QNMs can be found online at |http : //math . berkel eyT] 



edu/~dyatlov/qnmskds, We explain why the presented method gives the quantization 



conditions of Propositions L5 and |1.6[ but we do not provide a rigorous proof. 

We now compare the pseudopoles given by quantization conditions to QNMs for the 
Kerr metric]^] computed by the authors of [BeCaStj using Leaver's continued fraction 
method — see |BeCaSt[ Section 4.6] for an overview of the method and [BePaWi] Ap- 
pendix E] and |BeKo[ Section IV] for more details. The QNM data for the case of scalar 



10 



The results of the present paper do not apply to the Kerr case A = 0, due to lack of control on the 



scattering resolvent at the asymptotically flat spatial infinity. However, the resonances described by (0.2) 
are generated by trapping, which is located in a compact set; therefore, we can still make sense of the 
quantization condition and compute approximate QNMs. 
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0.2 
-0.285 1- 




[BeCaSt] 

-N = 2 



Figure 5. Comparison of order 2 approximation to QNMs with the data 
of [BeUaStj . Here I — 1, . . . , 4, k — —1,-1 + 1,1 — 1,1 (left to right), and 
m = (top) and 1 (bottom). 



perturbations, studied in this paper, computed using Leaver's method can be found online 
at |http : //www . phy . olemiss . edu/~berti/ qnms . html| 

Figure |5| com pares the second order approximation to QNMs (that is, solution to the 
equation (B.2) constructed using Qj and for j < 2) to the QNMs of [BeCaStj . Each 
branch on the picture shows the trajectory of the QNM with fixed parameters m, I, k for 
a G [0,0.25]; the marked points correspond to a = 0, 0.05, 0.25. The branches for 
same m, I and different k converge to the Schwarzschild QNMs as a — > 0. We see that the 
approximation gets better when / increases, but worse if one increases m\ this agrees well 
with the fact that the computed quantization conditions are expected to work when / is 
large and m is bounded. 

The left part of Figure [6] compares the second and fourth order approximations with 
the QNMs of [BeCa Stj (with the same values of a as before); we see that the fourth order 
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Figure 6. Left: comparison of order 2 and 4 approximations to QNMs with 
the data of [ BeCaSt] . Here I = 3, 4, I' — 0, 1, and m = 0. Right: log-log plot 
of the error of order 1-4 approximations to QNMs, as compared to [BeCaSt] . 
Here a = 0.1, k = I, m = 0, and I, plotted on the x axis, ranges from 1 to 7. 

approximation is considerably more accurate than the second order one, and the former is 
more accurate for a smaller value of I'. Finally, the right part of Figure [6] is a log-log plot 
of the error of approximations of degree 1 through 4, as a function of I; we see that the 
error decreases polynomially in I. 

B.2. Barrier-top resonances. Here we study a general spectral problem to which we will 
reduce both the radial and the angular problems in the next subsection. Our computation 
is based on the following observation: when the quantization condition of Section |4.3| is 



satisfied, the function u + has the microlocal form (4.3), with the symbol behaving like 
(r — ro) m near the trapped set. This can be seen from the proof of Proposition 4.6 if 
= —ihm, then u±(x) = x m and Biu± has to have the form (4.3). The calculations below 
are similar to |DiSj[ Section 3]. 

Consider the operator 

P v = D v A(y)D v + B(y;u,k). (B.4) 
Here the function A(y) is independent of u>, k, real-valued, and A(0) > 0; B is a symbol of 
order 2: 

B(y;u,k) ~ ^ Bj(y; u, k), 

j>0 
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with Bj homogeneous of degree 2 — j in the sense of (B.3). We also require that B be 
real-valued and 

B' (0;u,k) = 0, B£(0;w,fc) < 0. 

We will describe an algorithm to find the quantization condition for eigenvalues A of P y 
with eigenfunctions having the outgoing WKB form (B.6) near y = 0; we will compute A 
as a symbol of order 2: 

A ~ S2 k), Xj(M s u, sk) = s 2 ~ j X(u, k), s > 0. 
i>o 

More precisely, we will show how to inductively compute each Xj. The principal part Ao is 
given by the following barrier-top condition: 

A = B o {0;u,k). (B.5) 

In this case, we have 

B (y; u, k) = Ao(w, fc) - y 2 U (y; u, k), 

where U is a smooth function, and U (0) = — V"(0)/2 > 0. Define the phase function 
ip (y;u,k) such that 

^o(y; w, k) = yy/U (y;u; k)/A(y)\ 
note that ipo is homogeneous of degree 1. We will look for eigenfunctions of the WKB form 

u(y; u, k) = j^"^a{y\ u, k), (B.6) 

solving the equation P y u = Xu up to 0(\u\ + \k\)~°° error near y — 0. Here a is a symbol 
of order zero: 

a(y;u,k) ~ ^aj(y;uj } k) } 
with a,j homogeneous of order — j. 

Substituting (B.6) into the equation P y u = Xu and gathering terms with the same degree 
of homogeneity, we get the following system of transport equations: 

(L - B x + Ai)a 3 - = -LiOj-i + 2j (^+1 ~ \+i) a j-h j > °, 

0</<i 

L = 2i<4>' Ad y + z(A^)' = 2iy/U (y)A(y)yd y + i{yy/U {y)A{y))'; 

Li = d y A(y)d y , 

with the convention a_i = 0. 

Now, consider the space of infinite sequences 

C°° = {a = (a J )~ I a J G C} 
and the operator T : C°°(R) -)■ C°° defined by 

T(a) = a, &J - = aja(0)/i!. 



(B.7) 
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Let the operators L.,, B., : C°° — > C°° be defined by the relations 

TLj = L^X 1 , TBj = Bj-T 1 . 

We treat Lj, B,, as infinite dimensional matrices. We see that each Bj is lower triangular, 
with elements on the diagonal given by -Bj-(O); (Li)jk = for j + 2 < k. As for L , due to 
the factor y in front of the differentiation it is lower triangular and 

(L ) ij =i(2j + l)VC/ (0)A(0). 

One can show that there exists a smooth nonzero function do solving (Lq — Bi + Ai)ao = 
if and only if one of the diagonal elements of the matrix Lo — Bi + Ai is zero (the kernel of 
this matrix being spanned by Too)- Let m > be the index of this diagonal element; this 
will be a parameter of the quantization condition. We can now find 



Ai = Bi(0) - i(2m + l)y/U (0)A(0). (B.8) 

Now, there exists a nonzero functional f on C°°, such that f (a) depends only on a , ... , a m , 
and f vanishes on the image of L — Bx + Ai. Moreover, one can show that the equation 
{Lq — Bi + \i)a = b has a smooth solution a if and only if f{Tb) = 0. 

Take ao to be a nonzero element of the kernel of Lo — B\ + Ai; we normalize it so that 
f(Tao) = 0. Put a.j = Taj] then the transport equations become 

(L - Bi + Ai)aj = -Lia^-i + ^ (B /+ i - Aj+O^-j, j > 0. (B.9) 

0<l<j 

We normalize each aj so that f(aj) = for j > 0. The j-th transport equation has a 
solution if and only if the f kills the right-hand side, which makes it possible to find 

Aj+i = f ( - Ua.,-1 + B '+i a i-0' * > °- ( B - 10 ) 

0<Kj 



Using the equations (B.5), (B.8), (B.10), and (B.9), we can find all Aj and &j inductively. 



B.3. Radial and angular quantization conditions. We start with the radial quanti- 
zation condition. Consider the original radial operator 

P r . = D r {A r D r ) + V r {r; u, k), 
V r {r; u, k) = -A;\l + af{{r 2 + a 2 )co - akf . 



It has the form (B.4), with 

y = r-r , A{y) = A r , B{y; u, k) = V r {r; u, k). (B.ll) 
Here r is the point where V r achieves its maximal value, corresponding to the trapped 



point xq in Section 4.1 Now the previous subsection applies, with the use of the outgoing 
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microlocalization mentioned in the beginning of that subsection. Using (B.5) and (B.8) 
we can compute near a = 0, 

2ak(l- 9 AM 2 



r = 3M 



27M 2 
1 - 9AM 2 



1 - 



9M Rew 
2ak 



9M 2 Reu 



Z + ( a 2 (|A;| 2 + M 2 )), 
(Reoo) 2 + 0(a 2 (\k\ 2 + \oo\ 2 )), 



(B.12) 



Go + Q{ 



i(m + 1/2) + 



3V3M co 12 
V1-9AM 2 



0(1) for a = 0; 



reintroducing the semiclassical parameter, we get the formulas for T T in Proposition 1.5 



Now, we consider the angular problem. Without loss of generality, we assume that k > 0. 
After the change of variables y = cos 9, the operator Pe\v' takes the form 



P y = D y (l-y 2 )(l + ay 2 )D y 



(1 + a) 2 (au(l - y 2 ) - kf 



;i-y 2 )(i + m/ 2 ) • 

We are now interested in the bottom of the well asymptotics for the eigenvalues of P y , 



with the parameter /' from (B.l) playing the role of the quantization parameter m. The 



critical point for the principal symbol of the operator P y is (0,0). To reduce the bottom 
of the well problem to the barrier-top problem, we formally rescale in the complex plane, 
introducing the parameter y' = e t7T ^y, so that (y 1 ) 2 = iy 2 . We do not provide a rigorous 



justification for such an operation; we only note that the WKB solution of ( |B.6[ ) looks like 



e~ cy *a near y = for some positive constant c; therefore, it is exponentially 
decaying away from the origin, reminding one of the exponentially decaying Gaussians 
featured in the bottom of the well asymptotics (see for example |DiSj[ Section 3] or the 
discussion following [SaBaZwl Proposition 4.3]). There is a similar calculation of the bottom 
of the well resonances based on quantum Birkhoff normal form; see for example |CdVGui] . 
The rescaled operator P y i = —iP y takes the form (B.4), with y' taking the place of y and 

At >\ f-i , -r /\2\fi ■ i >\2\ ui * n i{l + a) 2 {au{l+i{y') 2 )-k) 2 

A(y ) = (1 + i(y ) )(1 - %a{y ) ), B(y;u,k) = (1 - - (y/)2)(1 _ - a(y/)2) • ( B -!3) 



We can now formally apply the results of Section B.2 note that, even though A and B are 
not real-valued, we have 

A(0) = 1, B (0) = -i(l + a) 2 (aReu - k) 2 , B' \0) < 0. 



An interesting note is that when a 

the spherical harmonics A = 1(1 + 1) exactly and without the assumption (B.l) 



and k > 0, the process described in Section |B . 2| gives 

In fact, 



the first three terms of the asymptotic expansion of A sum to 1(1 + 1) and the remaining 
terms are zero. 
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